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Abstract 


This  final  report  details  plasma-chemical  and  energy 
transfer  studies  carried  out  at  the  Research  Institute  for 
Engineering  Sciences,  Wayne  State  diversity.  Experimental 
and  analytical  research  has  been  aimed  at  understanding 
additive  effects  in  discharge  supported  molecular  laser  systems, 
particularly  the  CO  laser  system.  Extensive  numeric? 1 model 
codes  have  been  develooed  and  applied  to  the  CO  system  with 
considerable  success.  Experiments  in  our  laboratories  have 
provided  important  information  on  high  lying  vibrational  states 
in  CO,  in  agreement  with  model  predictions.  Collaborative 


studies  between  RIES  and  Princeton  University  on  the  behavior 


of  vibrational  distributions  under  strongly  pumped  conditions 


is  outlined. 
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Introduction 


This  report  summarizes  research  in  the  Research  Institute  for 

Engineering  Sciences  at  Wayne  State  University  into  enhanced  power 

and  efficiency  induced  by  additive  species  in  the  CO  molecular  laser 

system.  This  laser  system  is  characterized  by  high  efficiencies, 

60%  in  pulsed  mode^  and  up  to  50%  under  CW  conditions^,  significant 

(3) 

power  densities  in  a relatively  small  volume,  of  order  100  KW'  and 
a rather  broad  lasing  wavelength  range,  4.8  to  about  8 microns.  These 
features  are  often  a result  of  additive  species  which  significantly 
improve  laser  performance.  Many  reviews  of  the  CO  laser  system  have 
been  performed  in  recent  years including  our  own  unpublished 
review^.  In  this  final  report,  we  will  concentrate  on  the  two  major 
developments  from  our  research  program;  a detailed  study  of  the  effects 
of  additive  gases  on  laser  performance  and  collaborative  studies  with 
H.  Rabitz  and  S.  H.  Lam  at  Princeton  University  on  energy  transfer  processes 
in  the  upper  vibrational  states  of  the  CO  system.  The  majority  of  our 
work  to  date  in  these  areas  has  appeared  either  in  technical  report  or 
journal  form,  therefore,  this  final  report  will  serve  to  draw  together 
this  work  in  a coherent  framework. 

It  is  justified  to  say  that  the  general  behavior  of  the  CO  laser 
system  is  rather  well  understood.  Model  calculations  by  several  groups(8-13) 
have  been  able  to  adequately  reproduce  measured  performance  under  a wide 
variety  of  conditions.  However,  a significant  range  of  operating  con- 
ditions in  CO  laser  systems  involves  additive  gases  whose  concentrations 
are  quite  small  but  whose  effects  are  dominant  on  laser  performance.  For 

example,  the  high  reported  efficencies  for  CO  laser  performance  involving 
f 14 1 

Xe  and  Hg  additives'  . Although  the  role  of  these  gases  is  very  likely 
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to  enhance  ionization  characteristics  of  CO-He  and  promote  discharge 

stability,  other  additives  such  as  0^  and  are  not  amenable  to  such 

(151 

simple  interpretations  . 

In  general,  the  interest  in  our  studies  has  been  to  discharge  sup- 
ported low  pressure  CO  laser  systems.  Thus,  our  modeling  codes  involved 
the  early  development  of  a non-equilibruim  electron  energy  distribution 
code^).  This  code  has  been  instrumental  in  our  subsequent  energy  trans- 
fer and  plasma  chemistry  studies  as  it  has  provided  us  with  the  basic 
input  rates  for  modeling  and  simulation  studies.  It  is  worth  remarking  here 
that  in  molecular  or  atomic  systems  which  exhibit  strong  metastability  toward 
internal  energy,  i.e.  vibrational  and  in  which  large  energy  densities  are 

involved, superclastic  energy  feedback  from  the  lasing  species  to  the 

(151 

electrons  is  very  important  . Superelastic  processes  are  the  reverse  of 
electron  excitation  processes,  i.e.  reactions  in  which  the  electrons  gain 
energy  in  collisions  with  vibrationally  excited  CO  molecules.  These  processes 
are  important  in  discharge  environments  when  the  molecular  internal  energy 
states  (vibration)  are  highly  excited  as  in  the  CO  laser.  Thus,  electron 
processes  in  high  vibrational  states  can  be  important  as  an  energy  loss 
mechanism  even  though  the  reverse  excitation  process  is  small  compared  to 
other  energy  transfer  processes  in  these  states.  As  an  example  of  this  effect, 
e demonstration  calculation  on  the  shape  of  the  CO  vibrational  distribution 
both  with  and  without  upper  vibrational  state  electron  processes  is  shown  in 
Figure  1.  Conditions  typical  of  low  pressure  CO  discharge  laser  systems 
have  been  chosen  in  the  calculation.  The  details  of  the  numerical  output  for 
the  calculation  with  upper  state  processes  is  given  in  Table  I.  By  reference 
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Figure  1:  Comparison  of  a CO  vibrational, 

distribution  with  and  without  upper  state  e-processes. 
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Table  I:  Steady  State  Rates  in  CO 


(I-D 


STEADY  STATE  VIBRATIONAL  DISTRIBUTIONS  IN  CO 


TEMP  = 155. K TOTAL  PRESSURE  = 8.5  TOPR 


CC  DENSITY  = C 

.662E+1  7 

HE  DENSITY  = 

0.463F+-18 

O D 3NSTTY 

E DENSITY  = 0. 

200E* 1 1 

E/N  = 0 . 15P  E-1 5 

ELECTRON  TEMP 

= 13793. 

DRIFT  VELOCITY  = 0.371E+07 

V 

C (V) 

T (V) 

V-V* 

GAIN 

JM  A X 

0 

4. 648E+16 

2500. 

0-  0 

0.0 

0 

1 

1. 135E+16 

2189. 

1-  0 

0.0 

0 

2 

3 . 54 1 E* 15 

26  17. 

2-  1 

7. 214E-08 

29 

3 

1 . 406E+15 

32  58. 

3-  2 

6.323  £-"6 

24 

4 

7. 028  £♦ 1 4 

4288. 

4-  3 

2. 085E-04 

19 

5 

4.  308E  + 1 4 

6995. 

5-  4 

2 . 2 84  E-r  3 

15 

6 

3.092E+14 

8730  . 

6-  5 

9 . 2C8  E-°  3 

1 1 

7 

2.  46 1 E* 1 4 

12837. 

7-  6 

1 .801E-02 

10 

8 

2.079E+14 

16741. 

8-  7 

2.405E-02 

9 

9 

1.  8 1 2 E*  1 4 

20231  . 

9-  ft 

2. 750E-02 

8 

10 

1.602E+1 4 

22274. 

10-  9 

2.83RE-02 

8 

11 

1. 424E+14 

23057. 

11-10 

2.790E-02 

8 

12 

1 . 268E+14 

23"00. 

12-11 

2.f 63E-02 

8 

13 

1.  128E  + 14 

22501  . 

13-12 

2.493  E-02 

8 

14 

1. 002E+14 

21953. 

14-13 

2.338E-02 

8 

15 

8. 885E+13 

21319. 

15-14 

2.146  E-02 

8 

18 

7. 870E+13 

20812. 

16-15 

1.963  E-02 

8 

17 

6. 963E+13 

20336. 

17-16 

1 . 792E-02 

9 

18 

6. 156E+13 

19909. 

18-17 

1 .634E-02 

9 

19 

5. 442E+13 

19603. 

19-18 

1 .4  89  E-0 2 

9 

20 

4 . 805  E* 1 3 

19108. 

20-19 

1.  3 26 E-0 2 

9 

21 

4. 238E+13 

18657. 

21-20 

1 .177E-02 

9 

22 

3. 730E+13 

18082. 

22-21 

1 .034E-02 

9 

23 

3.  279 E + 1 3 

17653. 

2 3-22 

9 .083 E-0 3 

9 

24 

2. 880E+13 

17223. 

24-23 

7. 941 E-0 3 

9 

25 

2. 523E+13 

16653  . 

25-24 

6. 8 57 E-0 3 

9 

26 

2. 206E+1 3 

16114. 

26-25 

5.890E-C  3 

9 

27 

1. 922E+1 3 

15474. 

27-26 

4 .991  E-0 3 

9 

28 

1.673E+13 

15071. 

28*27 

4 . 277  E-0  3 

10 

29 

1. 450E+13 

1 4437. 

29-28 

3 ,6rt4E-03 

10 

30 

1 . 253E+1 3 

13875. 

30-29 

2. 8 78  E-0 3 

10 

31 

1.077E+13 

13179. 

31-30 

2. 3 64  E-0  3 

10 

32 

9.  1ft  2 E ♦ 1 2 

12224. 

32-31 

1 .869E-03 

10 

33 

7. 811E+12 

11883. 

33-32 

1.542E-03 

10 

34 

6.623E+12 

11435. 

34-33 

1. 256E-03 

11 

35 

5. 596E+12 

10983. 

35-34 

1.022E-03 

1 1 

36 

4. 714E+12 

10590. 

36-35 

8. 3C7E-04 

1 1 

37 

3.  986  E*1 2 

10619. 

37-36 

7.^71 E-04 

1 1 

3ft 

3 . 265E* 1 2 

8758. 

38-37 

4. 524E-04 

12 

39 

2. 709E>  1 2 

9166. 

39-38 

4 . r 21 E-04 

12 

0.  1 00  E*C  6 
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ELECTRON  PATES 


FORM  ATICN 

PROCESSES 

LOSS  PROCESSES 

V 

(V-L)  TO  V 

( V*L)  TO  V 

V TO  (V-L) 

V TO  (V ♦ L) 

NET  PATE 

0 

0.  0 

0. 246E*1 9 

0.0 

0.856E+1 9 

-0.6  10E+19 

1 

0. 557F+19 

0.  1 6 5 E* 1 9 

0. 206F+19 

0 . 41 8E  + 1 9 

0.978E+18 

2 

0.416E+19 

0. 10RF+1 9 

C.  156E*  1 9 

0. 196E  + 1 9 

0. 172E+19 

3 

0. 265E+19 

0. 803E+1  8 

0. 106F+19 

0 . 104E+1 9 

0.  136E+19 

u 

0.  173E+19 

0.677E+18 

0. 776E+18 

0.647E+18 

0.  983E*  1 8 

5 

0.  1 16E+19 

0. 624E* 1 8 

0.639F+18 

0.476E+18 

0. 664E+18 

€ 

0. 859E+18 

0. 603E+18 

0. 584E+18 

0.399F+18 

0.480E+18 

7 

0. 650E+18 

0. 595E+1 8 

0.568E+  18 

0.363E+1 8 

0.  3 14E+18 

fl 

0. 498E+ 18 

0. 588E+18 

0. 56  9F*  18 

0. 345E+18 

0. 172E+ 1 8 

9 

0. 405E+18 

0. 579E+1 8 

0.573Et18 

0.334E4-1  8 

0.779E+17 

10 

0.  36 2E+- 1 8 

C.566E+18 

0. 574E+  18 

0.325E+18 

0. 292E+ 17 

11 

0. 339E+18 

0.  54  9 E* 1 8 

0.571E+18 

0.315E4-1  8 

0. 257 E* 16 

12 

0.325E+18 

0. 529E+ 1 8 

0.563F+18 

0. 304E  + 1R 

-0. 131E+17 

13 

0. 312E+18 

0.  505E+18 

0 . 54  8E  ♦ 1 8 

0.291E+18 

-0.  223E+17 

1U 

0. 3C0E+18 

C.479E+18 

0. 530E+  18 

C.277E+18 

- 0.  2 80  F+  1 7 

15 

0. 286E+ 18 

0. 452E«-1  8 

0.508E+18 

0.262E+18 

-0.3 11E+17 

16 

0. 272E+18 

0. 425E+1 8 

0. 483E+- 18 

0.246E+18 

-0. 328E+17 

17 

0. 257S+18 

0. 397E+1 8 

0.457E+18 

0.231E+1 8 

-0. 336E+17 

18 

0.  242E-*- 1 8 

0.  370F+18 

0. 430E*18 

0.215E+1 8 

-0.  338E+  17 

19 

0.  227E+18 

0.  34 3 £■»•  1 8 

0.404E+18 

0 .20 1F+18 

-0. 341  E*  17 

20 

0.212E+18 

0. 317E+18 

0. 377E+18 

0. 186E  + 18 

-0.  337E+17 

21 

0.  197E+18 

0.  292F.  + 1 8 

0.350E+18 

0.172E+18 

-0.  3 3 2 E+  1 7 

22 

0.  1 82E*  1 8 

C.268E+18 

0. 324E+18 

0. 158E+18 

- 0.  320E+17 

23 

0.  168E+18 

0. 245E+1 8 

0.299E+18 

0.145E+18 

-0.  308 E*  17 

24 

0.  155S+18 

C. 223E+18 

0. 275E+18 

0. 133E  + 18 

-0.298E+17 

25 

0.  142E+18 

C.  202E  + 18 

0.252E+  18 

0.121E+1 8 

-0.  284E+17 

26 

0.  130E+18 

0. 182E+18 

0.  229E+  18 

0. 1 10E  + 18 

-0.  27 1 E+  1 7 

27 

0.  118E+18 

C. 163E+1 8 

0.208E+18 

0.991 E + 1 7 

-0.  253E4-17 

28 

0.  107E* 18 

0. 146F+1 8 

0.  188E+18 

0.894E+17 

-0.241E+17 

29 

0.  97  IE*  1 7 

0.  130E+18 

0. 169E  + 18 

0.801E+17 

-0.  2 2 5 E* 1 7 

30 

0.874E+17 

0. 115E+18 

0. 152E-HB 

0.716E+17 

-0. 211E+17 

31 

0.  784E-M7 

0. 101E+18 

0. 135E+18 

0 .635E  + 1 7 

-0. 196E+17 

32 

0.699E+17 

0.  875  £♦  1 7 

0. 119E+18 

0. 557E+17 

-0.  172E+17 

33 

0.618E+17 

0.748  F.  *17 

0. 104E+18 

O .U83E  + 1 7 

-0.  161E+17 

34 

0.545E+17 

0.629E+17 

0. 914F*17 

0.414E+17 

-0. 154E+17 

35 

0.  478E+17 

0.521E*17 

0. 796E+17 

C.350E+17 

-0. 148E+17 

36 

0.417E+17 

0.414E+17 

0.691F+17 

0.288E+17 

-0.  1 47F>  1 7 

37 

0.  363E+17 

0. 305E+17 

0.601E+17 

0.228E+1 7 

-0.  16 2 E* 1 7 

38 

0.  315E+17 

C.  192E+1 7 

0.  506E*  17 

0. 153E  + 1 7 

-0.  152E+17 

39 

0.  269E+17 

0.0 

0.431F+-17 

0.0 

-0. 162E+17 

NET  ELECTRON  ENERGY  INTO  CO  = 4.871E+22 

NET  SUPERELASTIC  ENERGY  OUT  OF  CO  = -2.  960E*22 
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(1-3) 


VIERATICN-VIBRATION  PATES 


FORMATION 

PPO CESSES 

loss  : 

PROCESSES 

V 

(V-1)  TO  V 

(V+1)  TO  V 

V TO  (V-1) 

V TO  (V  + 1) 

NET  RATE 

0 

0.0 

0.747E+21 

n.  n 

0.7422+21 

C.  5 7 0 E+  1 9 

1 

0. 729E+21 

0. 742E+21 

0. 747E+21 

0.7242+21 

-0.  8 392+  18 

2 

0. 400E+21 

0. 724E+21 

0. 729P+21 

0.3972  + 2 1 

-0.  1632+19 

3 

0.  178E+21 

0. 397E+21 

C. 4002+21 

0. 1762+21 

-0.1312+19 

4 

0. 817E+20 

0. 176E+21 

0. 178E+21 

0.808E+20 

-0. 96  8 E+  18 

5 

0. 449E+20 

0. 808E+20 

0. 817E+20 

0.446E+20 

-0.  660E+18 

6 

0. 301E+20 

0. 446E+2 0 

0.44  9F  + 2 0 

0. 3032  + 20 

-0. 4882+18 

7 

0.  232E+20 

C. 303E+2P 

C.  301E+20 

0.237^+20 

-0.  3222+  1 8 

8 

0.  197E+2C 

0. 237E+20 

0. 232E+20 

0.204E+20 

-0.  189E+  1 8 

9 

0.  175E+20 

0. 204E+20 

0.  197E  + 20 

0. 1832  + 20 

- 0.  9502+17 

10 

0.  159E+2C 

C.  183E  + 20 

0. 175E+20 

O.168E+20 

-0.488E+17 

11 

0.  146E+20 

0.  168E  + 20 

0. 159E+20 

0.  1 5 CE  +2  0 

-0.  1702+  17 

12 

0. 134E+20 

0.  155E  + 20 

0. 146E+20 

0 . 1 432  + 20 

-0.692E+16 

13 

0.  1232*20 

0. 143E+20 

0. 134E+20 

0. 132E+20 

-0.  1 36E+  16 

14 

0. 1 12E+20 

0.  132E  + 20 

0. 123E+20 

0.121E+20 

0.  8192+16 

15 

0. 1C2E+20 

0. 121E+20 

0.11 2E  + 20 

0.11 1E+20 

0. 156E+17 

16 

0. 925E+19 

0.  1 1 1E  + 20 

0. 1022+20 

0 . 1 0 1E  + 2 0 

0.  226E+17 

17 

0.  837E+1  9 

0.  101E  + 20 

0. 925E+19 

C .9172+1 9 

0.  2802+  17 

18 

C. 756E+19 

0. 917E+1 9 

0.837E+19 

0 . 8 3 4 E ♦ 1 9 

0.  3202+17 

19 

0. 6812+19 

0. 834E+1 9 

0.756E+19 

0.7662+19 

0.  258E+  17 

20 

0.61’  19 

0. 756E+1 9 

0. 6812  + 19 

0.684E  + 1 9 

0.  3032+17 

21 

C.  c 9 

0. 6P4E+19 

0. 613E+19 

0.61 9E +19 

0.  2602+17 

22 

0 19 

0. 61  9F.+  1 9 

0.550E+19 

0.5582+1 9 

0.  3 35E+17 

23 

0. 558E+19 

0. 493E+19 

0.50  32  + 19 

0.3242+17 

24 

+ 19 

0. 503E+19 

0.441E+19 

0.4532+19 

0. 2652+17 

25 

j.  >4  9S+ 1 9 

0. 453E+1 9 

0.393E+19 

0.406E+1 9 

0.  2762+17 

26 

0. 3C9E+19 

0. 406E+19 

0. 3492+19 

0. 3632+19 

0. 243E+17 

2 r 

0. 272E+19 

0. 363E+1 9 

0. 309E+19 

0.323E+19 

0.  3 17E+17 

28 

0. 237E+19 

C. 323E+1 9 

0.2722+19 

0.2862+1 9 

0. 236E+ 17 

29 

0. 2C6E+19 

0. 286E+19 

0.237E+19 

0.2532+19 

0.  247E+  17 

30 

0.  177E+19 

0. 253E+19 

0. 2062+19 

0.222E+19 

0.  201E+17 

31 

0.  150E  + 19 

0. 222E+1 9 

0.  1772+1  9 

0.193E+19 

0.1432+17 

32 

0. 125E+19 

0.  193E+19 

0. 150E+19 

0. 166E  + 1 9 

0. 277E+ 17 

33 

O'.  103E  + 19 

0.  166E  + 1 9 

0.  12SF+  19 

0.1412+19 

0. 215E+17 

34 

0. 828E+ 18 

0.  141E+19 

0.  1032+  19 

0. 1 19E+19 

0.  184P+17 

35 

0. 645E+1 8 

0. 1 19E+1 9 

0.828E+18 

0.9912+18 

0. 173 E+ 17 

36 

0.480E+18 

0.991E+18 

0. 6452+18 

0.8052+1 8 

0.2  17E+17 

37 

0. 396E+18 

0. 8C5E+1 8 

0.4802+18 

0.7002+1 8 

0. 2132+17 

38 

0.296E+18 

0.7’0E+18 

0. 3962+18 

0.6R7E+18 

0. 125E+17 

39" 

0. 5872+18 

0.0 

0.296F+18 

C.O 
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the  second  page  of  Table  I,  labeled  "electron  rates"  we  note  that  in  vibra- 
tional states  above  v = ll  the  net  electron  rates  are  negative  indicating  a 
net  flow  of  eiergy  from  the  vibrational  levels  to  the  electrons.  This  feature 
is  characteristic  of  all  of  our  distribution  calculations  although  the 
quantitative  rates  depends  on  the  particular  cross  sections  used  in  the  upper 
states.  Also  in  the  CO  system  our  studies  have  indicated  that  electronic 
excitation  by  electrons  is  important  both  for  plasma  heating  as  a result  of 
subsequent  quenching  collisions  but  also  for  carbon  formation.  Although 
carbon  is  not  predicted  to  be  important  in  the  plasma  chemistry  of  the  CO 
discharge  laseP5^,  it  does  limit  the  operating  time  due  to  deposits  which 
increase  cavity  losses  leading  eventually  to  a reduction  in  lasing  power 
output.  Therefore,  the  electronic  cross  sections  which  are  rather  poorly  known 

(15) 

become  important  in  laroe  scale  CO  predictions  '. 

In  modeling  the  dynamic  and  steady  state  behavior  of  the  CO  system 
we  have  also  developed  and  applied  time  dependent  and  C W codes  incorpor- 
ating stimulated  emission  and  state-of-the-art  rate  coefficients.  These 
codes  have  enabled  our  understanding  of  CO  laser  performance  to  become 
quite  sophisticated  in  those  cases  where  the  basic  molecular  rate  information 
is  available.  In  addition,  these  model  calculations  have  indicated  priorities 

on  important  rate  processes  for  which  poor  information  is  currently  avai Table. 
Our  studies  on  0,  effects  in  the  CO  laser  plasma  and  subsequent  effects  on 
laser  performance  as  detailed  in  Appendix  I have  indicated  the  importance 
of  plasma  chemistry  effects  in  low  pressure  molecular  discharges  and  rein- 
forced the  need  for  continued  work  in  providing  basic  rate  information  on 
molecular  systems  as  a precurson  to  our  understanding  and  hence  optimization 
of  laser  systems. 
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Our  understanding  of  the  detailed  energy  transfer  processes  in  the 
CO  system  was  given  a significant  contribution  as  a result  of  the  vibration- 
translation  calculations  on  CO-He  by  Ver ter  and  Rabitz' Collaboration 
between  the  groups  in  RIES  and  Princeton  led  to  a series  of  papers  at  the 
5th  Conference  on  Chemical  and  Molecular  lasers  in  St.  Louis  in  April 
(1977)-see  Appendix  II.  The  collaboration  was  particularly  successful  as 
it  blended  the  theoretical  and  analytical  effort  at  Princeton  with  the 
experimental  and  modeling  capabilities  in  RIES. 

It  has  been  known  for  a long  time  that  energy  transfer  processes  in 
the  upper  vibrational  states  of  CO  play  an  important  role  in  CO  laser  be- 
havior. Substantial  output  power  is  emitted  from  vibrational  states  between 

v-10  and  v-15  while  even  higher  vibrational  states  are  important  in  the 
(15) 

dissociation  of  CO  . Yet  until  recently,  no  direct  calculation  or 
experimental  measurements  have  been  available  to  incorporate  into  laser  codes 
for  performance  predictions.  The  Verter-Rabi tz  calculations  provided  for  the 
first  time  some  sound  calculations  on  the  temperature  and  level  dependence 
of  V T rates  in  CO-He.  It  is  significant  to  note  that  large  deviations  were 
noted  for  upper  vibrational  levels,  i.e.  above  about  v=30,  from  the  previously 
accepted  SSH  results^’^.  These  conclusions  are  detailed  in  Appendix  III^) 
The  experimental  technique  developed  during  this  program  involves 
direct  observation  of  the  spontaneous  end-light  emission  on  overtone  tran- 

( 0 

sitions  in  a CO  laser  plasma.  The  original  work  within  RIES  on  this  technique 
is  given  in  Appendix  IV.  Currently,  we  are  able  to  observe  and  interpret 
the  CO  vibrational  distribution  up  to  v=35  under  varying  plasma  chemical 
and  physical  conditons.  The  observation  of  the  entire  vibrational  dis- 
tribution through  the  second  "knee"  provides  a unique  comparison  to  both 
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(18)  (19) 

model'  ' and  analytical  ' interpretations  of  energy  transfer  in  high 
vibrational  levels.  These  collaborative  studies  are  continuing  to  provide 
insight  into  other  gase  mixtures,  and  a wider  range  of  temperatures. 

Lastly,  our  studies  have  pointed  out  the  importance  of  plasma-chemical 
processes  in  discharge  supported  CO  lasers'  ' . These  extensive  model  cal- 
culations suggested  that  rather  small  amoun  ts  of  0-atoms,  particularly 
produced  by  direct  dissociation  of  0^  by  electrons,  can  have  a detri- 
mental effect  on  laser  performance  due  to  upoer  state  V-T  processes. 

Recent  experimental  measurements  in  our  laboratory  on  vibrational 
distributions  in  a CO  laser  in  the  presence  of  0^  have  collaborated  these 
model  calculations''"0).  The  report  discribing  these  measurements  and 
detailing  the  conclusions  can  be  formed  in  Appendix  V. 
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ABSTRACT 

Trace  additions  of  several  gases  are  known  to 
significantly  effect  CO  laser  power  and  efficiency. 

Since  small  concentrations  of  additives  are  generally 
involved,  these  additives  are  recognized  as  predomin- 
antly influencing  the  characteristics  of  the  laser 
plasma.  Oxygen  has  been  an  additive  receivino  substantial 
attention  in  past  measurements.  In  this  paper  we  develop 
a consistent  interpretation  for  all  the  observed  oxygen 
effects  both  from  our  laboratories  and  elsewhere.  Oxygen 
in  small  concentrations  will  be  shovm  to  improve  laser 

- - 

efficiency  and  power  through  its  effect  on  the  electron- 
ion  recombinations  process.  At  higher  oxygen  additions, 
laser  output  is  degraded  due  to  vibration-translation 
relaxation  by  oxygen  atoms.  In  addition,  carbon  formation 
and  removal,  plasma  heating  and  cooling  effects,  hetero- 
geneous loss  processes  and  the  influence  of  polymer  ions 
will  be  discussed  through  model  calculations. 


| Oxygon  as  an  additive  in  discharge  sustained  CO 


molecular  laser  systems  has  been  used  extensively  to 
improve  laser  output.  Although  numerous  observations  on  the 
role  of  oxygen  have  been  made,  no  consistent  interpretation 
has  been  developed  which  characterizes  these  observations. 

In  this  paper  we  review  the  observations  on  oxygen  as  an  ad- 
ditive and  present  a detailed  model  which  provides  a consis- 
tent framework  for  interpreting  these  observations  as  well  as 
the  effects  of  other  additives  in  the  CO  laser  system. 

In  order  to  accomplish  this  understanding,  the  character- 
istics of  the  electron  energy  distribution  have  been  studied 
through  a numerical  solution  to  the  relevant  Boltzmann  equa- 
tion and  the  detailed  plasma  and  neutral  chemistry  and  energy 
transfer  processes  have  been  characterized  through  a time  de- 
pendent rate  equation  model.  These  models  are  detailed  in 
this  paper  with  application  to  the  CO/02/He  system.  Two  major 
points  are  indicated  as  a result  of  this  study  on  the  effect 
of  oxygen  on  CO  laser  performance:  ionic  and  neutral  chemistry 
are  very  important  to  an  understanding  of  the  observations  on 
oxygen  effects,  and  superelastic  processes  are  critical  to 
explaining  observed  ionization  levels. 
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A uniform  feature  of  additive  gases  found  to  enhance 
laser  power  and  efficiency  is  the  fact  that  these  gases  have 
lower  ionization  potentials  than  either  CO  or  He.  Thus,  the 
additive  can  generally  be  viewed  as  being  a source  of  electrons 
at  an  E/N  (electric  field/total  number  density)  more  appropriate 
to  maintaining  the  mean  electron  energy  near  that  for  optimum 
coupling  to  the  vibrational  states  of  CO.  Although  this 
viewpoint  is  undoubtedly  correct  for  some  additives  of  general 
interest,  e.g.  Xe , oxygen  will  be  shown  to  enhance  laser  output 
and  efficiency  as  a result  of  its  effect  on  the  ion  recombination 
dynamics  in  the  plasma. 

In  this  paner,  homogeneous  plasma  chemistry  and  energy 
transfer  processes  will  be  stressed.  At  total  pressures  around 
10  torr  under  liquid  nitrogen  cooling,  it  will  be  shown  that 
ambipolar  diffusion  competes  about  eaually  with  homogeneous 
recombination  as  the  major  electron  loss  process.  At  higher 
pressures,  the  electron  density  is  dominated  by  ionization  and 
homogeneous  recombination.  Attachment  processes  will  be  shown 
to  be  generally  unimportant  due  to  the  large  detachment  rate 
by  CO. 

II . ELECTRON  ENERGY  DISTRIBUTION  CALCULATIONS 

A.  Electron  Impact  Cross  Section  Data. 

A summary  of  the  electron  impact  cross  section  data  used 
in  this  study  appears  in  Table  1.  Special  mention  should  be 
made  of  the  cross  sections  for  electronic  excitation  of  CO. 

Near  threshold  there  is  little  information  on  the  magnitude  of 
the  cross  section  for  electronic  excitation  in  CO.  Most  pre- 
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TABLE  1 

ELECTRON  IMPACT  PROCESSES  IN  THE  C0-02  SYSTEM 


Electron  Impact  Processes  Energy  Loss  (eV)  Reference 


e + 

CO: 

1. 

CO^  (8  levels) 

0.266  - 2.034 

2 

2. 

CO(a3n) 

6.04 

see 

text 

3. 

CO(A1n) 

8.07 

see 

text 

4. 

c + o" 

9.00 

3 

5. 

C +0 

12.55 

3 

6. 

CO+ 

14.013 

4 

e + 

—2  : 

7. 

0+  (8  levels) 

0.193  - 1.46 

4 

8. 

0-  (a2A  ) 
2 g 

0.98 

5 

9. 

°9  <blO 
2 g 

1.64 

5 

10. 

0_ (A3£+) 

4.5 

1, 

6 

2 u 

11. 

°2(b3eu> 

8.4 

1, 

6 

12. 

9.7  eV  allowed 

9.7 

6 

13. 

o + o“ 

3.62 

3 

14. 

12.063 

4 

15. 

o 

+ 

+ 

o 

19.54 

3 

e + 

_0 : 

16. 

0(1D) 

1.96 

7 

17. 

0(1S) 

4.17 

7 

18. 

0+ 

13.6 

4 

_ 

e + 

°9: 

— Z 

19. 

°2 

11.0 

see 

text 

Note  that  t refers  to  vibrational  excitation  and  * refers 
tc  electronic  excitation. 
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vious  workers  have  used  a single  composite  cross  section,  due 
to  Hake  and  Phelps/  ' with  a threshold  of  6 eV  and  a peak  vul 
of  5 x lO"'*'6  cm2  at  10  eV  corresponding  to  the  C0(a2Il)  state. 

This  cross  section,  however,  is  subject  to  considerable 
uncertainty  as  pointed  out  by  Hake  and  Phelps. 


For  total  cross  sections  from  threshold  up  to  about  25 

( 8 ) 3 

eV,  the  work  of  Ajello  on  the  excitation  of  CO  (a  II)  and 

CO(A^II)  has  been  used  together  with  the  theoretical  work  of 

Chung  and  Lin  J 'on  the  excitation  of  eleven  CO  electronic  states. 

In  the  calculations  reported  here  we  have  included  the  two  CO 

electronic  states  with  the  largest  cross  sections:  the  a2n 

state  at  6.04  eV  and  the  A ^ IT  state  at  8.07  eV.  The  cross 

sections  used  are  composites  of  those  given  in  References  8 

and  9 . They  arc  in  substantial  agreement  with  a similar  set 

(10) 

of  cross  sections  due  to  Sawada,  el;  al.  y 

3 3 

With  reference  to  the  exciting  of  the  A I and  B E 
states  of  O2 , it  is  known  that  these  excited  states  pre- 
dissociate. The  (^(A^E^)  state  predissociates  into 

0(^PT  + 0(^P)  and  92(B^Eu)  predissociates  into  0('^P)  + 

0 (1D)  . ^1:L  ’ 12-)  Although  there  are  several  ways  for  the 

3 3 

02(B)  state  to  ^redissociate  into  0(  P)  + 0(  P)  via  curve 

(12) 

crossings,  the  probabilities  of  this  occurring  are 

unknown.  It  has  been  assumed  in  this  work  that  all  of  the 

3 1 

excitation  follows  the  P + D predissociative  path.  The 
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•k 

cross  section  for  ionization  of  the  O2  metastable  is 
unknown  and  has  been  taken  to  be  the  same  as  that  for 
ionization  of  O2 (X)  displaced  in  energy  by  the  excitation 

k 

energy  of  . This  is  about  the  least  arbitrary  method 

of  obtaining  an  estimate  of  the  cross  section  for  this 

process  in  the  absence  of  experimental  or  theoretical  data. 

There  is  some  support  for  this  assumption  from  the  work  of 

Burrowv  ' on  dissociative  attachment  from  . He  measured 

the  ratio  of  the  ionization  cross  section  of  the  a^A  state 

g 

3 

to  that  of  the  X state  at  approximately  0.5  eV  above 
each  respective  threshold  and  found  it  to  be  0.8  ± 0.3. 

B.  Results  for  CO. 

Nighan^^  has  shown  through  electron  energy  distribution 
calculations  that  as  the  E/N  of  a discharge  is  increased,  the 
fraction  of  energy  going  into  electronic  excitation  of  CO 
increases  at  the  expense  of  vibrational  excitation  suggesting 
that  laser  operation  at  low  E/N  is  desirable.  This  is  il- 
lustrated in  Figure  1 which  plots  the  fraction  of  electron 
energy  going  into  the  CO  vibrational  and  electronics  states 
for  two  values  of  E/N.  Here  u^  = 2/3u,  where  u is  the  mean 
electron  energy.  The  effect  of  superelastic  collisions  with 
vibrationally  excited  CO  is  also  shown  in  the  Figure  1.  Elec- 
trons in  the  1-4  eV  energy  range  gain  quanta  of  energy  from 
the  excited  vibrational  states,  ranging  from  .27  eV  for  v = 1 
to  2.30  eV  for  v = 8,  increasing  the  mean  energy  of  the  dis- 
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FIGURE 


VIBATIONAL  LEVEL  LEVEL 


1:  Fractional  energy  balance  in  pure  CO  as  a function 

of  E/N.  The  effect  of  including  superelastic  processes 
is  demonstrated  by  the  circled  points. 

f 

* 
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tribution.  This  has  only  a slight  effect  on  the  vibrational 
excitation  rate  coefficients  but  dominates  the  electronic 
excitation  and  ionization  rate  coefficients  due  to  the  signi- 
ficant increase  in  the  number  of  electrons  in  the  high  energy 
tail  of  the  distribution  function.  This  superelastic  feed- 
back will  be  shown  later  to  be  of  great  importance  in  CO 
laser  operation. 


Ill . THE  CO-Op-He  DISCHARGE 


A Review  of  Experimental  Work  on  the  Effects  of  on  the 
CO  Laser. 

A number  of  experimenters^^  ^^have  found  that  a small 
amount  of  added  to  a CW  CO  laser  plasma  has  a pronounced 
effect  upon  the  operation  of  the  laser.  Bhaumik,  et.  alV  have 
found  that  small  amounts  of  O2  (approximately  5%  of  the  CO  par- 
tial pressure)  enhanced  the  power  output  by  a factor  of  10-20% 
and  reduced  carbon  deposits  on  the  walls  of  the  laser  tube. 

They  suggested  that  the  O^  affected  the  dissociation  reaction 
CO  t C + O by  driving  it  to  the  left.  Larger  amounts  of  O^ 
had  a deleterious  effect  on  laser  operation,  possibly  due  to 
the  large  electron  attachment  cross  section  of  O^  or  to  the 
formation  of  CO^,  in  their  opinion.  Hartwick  and  Walder  ' 
found  that  the  total  laser  output  increased  and  that  the  laser 
became  more  stable  with  the  addition  of  C>2  to  the  CO-He  dis- 
charge. The  optimum  mixture  was  found  to  be  28  torr  He, 

2 torr  CO  and  0.1  torr  02-  In  addition,  they  measured  the 
gas  temperature  in  the  plasma  using  a thermocouple  and  cor- 
related it  to  laser  output  finding  that  the  minimum  in  tem- 
perature, approximately  295°K,  corresponded  to  the  maximum  in 
power  output  as  the  O ^ flow  rate  was  varied  for  a liquid  ni- 
trogen cooled  laser.  At  slow  oxygen  flow  rates  the  measured 
kinetic  temperature  was  approximately  320°K. 
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Extensive  work  involving  the  effects  of  0 2 upon  CO  lasers 

(13  191 

has  been  done  by  Keren,  Avivi  and  Dothan.  ’ ' They  found 

that  the  addition  of  0^  increased  the  power  output  and  de- 
creased the  operating  E/N  of  a water  cooled  CO  flow  laser. 

At  a constant  current  of  17  mA  the  laser  output  power  reached 
a maximum  of  4 watts  with  the  addition  of  50  m-torr  of  02  to 
a mixture  of  1.4  torr  CO  in  18  torr  He  and  the  axial  field 
strength  in  the  discharge  dropped  from  136  V/cm  to  88  V/cm. 
Power  output  then  decreased  with  the  addition  of  larger  con- 
centrations of  O ^ and  lasing  was  completely  quenched  at  an  02 
partial  pressure  of  100  m-torr.  They  also  pointed  out  that 
the  discharge  was  very  unstable  and  that  laser  action  could 
not  be  obtained  in  the  absence  of  the  additive.  On  the 
basis  of  the  measured  decrease  in  E/N  and  a calculation  of  the 
fraction  of  electron  energy  flowing  into  CO  vibrational  and 
electronic  states  as  a function  of  E/N  (similar  to  that  pre- 
sented in  Figure  1 here)  they  concluded  that  the  enhanced 
power  output  was  due  to  more  efficient  vibrational  pumping  be- 
cause of  the  lower  E/N. 

Keren,  et  a_l^^also  analyzed  the  composition  of  the  ions 
leaving  the  discharge  using  a mass  spectrometer  and  found  that 
with  the  addition  of  20  m-torr  02  to  the  laser  plasma  (1.4% 
of  the  CO  partial  pressure) , 02  became  the  dominant  positive 
ion  in  the  discharge.  With  no  oxygen  present  the  discharge  was 
dominated  by  the  dimer  ion  and  higher  order  polymer  ions. 
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They  suggested  that  the  dominance  of  0^+  was  due  to  ionization 
via  a three  step  process: 


e~  + CO  -*•  CO*  (a J ) + e“ 

CO*  + 02  -*•  CO  + 0*  (metastable) 

0*  + e •>  0 2 + 2e~ 

claiming  that  the  probability  of  producing  0.,+  by  direct  ioni- 
zation is  at  least  a factor  of  twenty  lower  than  the  probabil- 
ity of  the  above  process.  They  dismiss  the  charge  transfer  re- 
action, 

C0+  + o2  -*•  o2+  + CO 

as  a means  of  forming  02+  because,  in  their  opinion,  it  cannot 
explain  the  decrease  in  E/N  that  occurs  with  the  addition  of 
molecular  oxygen.  Using  this  three  step  model,  they  compute 
that  about  5%  of  the  02  molecules  are  excited  to  0*  for  a 
total  02  partial  pressure  of  25  m-torr. 

Because  very  small  amounts  of  molecular  oxygen  added  to 
a CO  laser  have  significant  effects  upon  laser  performance, 
it  appears  clear  that  the  plasma  properties  of  the  discharge 
are  affected  but,  until  now,  a consistent  interpretation  has 
not  been  available.  The  calculations  that  follow  are  directed 
toward  providing  this  systematic  explanation  of  the  observed 
properties  of  additive  02  in  a CO  discharge. 
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B.  Plasma  Chemistry  Calculations 


In  order  to  assess  the  importance  of  some  of  the  processes 

that  can  occur  in  the  discharge,  the  computer  code  for  calcu- 

(20) 

lating  the  electron  distribution  function  was  combined  wita  a 
time-dependent  chemical  kinetics  code.  This  allows  the  con- 
centrations of  all  important  chemical  species  to  be  computed 
as  a function  of  time.  Since  we  are  studying  CW  laser  behav- 
ior, this  computation  is  carried  out  for  times  of  the  order  of 
contact  times  in  flowing  CW  laser  systems  (i.e.,  tenths  of 
seconds) . On  these  time  scales  many  of  the  species  concentra- 
tions are  in  steady  state.  The  differential  rate  equations  are 
integrated  directly  rather  than  solving  the  steady  state  equa- 
tions because  the  coupled  non-linear  algebraic  equations  that 
arise  from  a steady  state  analysis  are  difficult  to  solve  in 
a systematic  way.  Uore  importantly , there  are  some  species  in  a 
.fast  flow  laser  (such  as  C02  in  a CO  laser)  that  may  not  reach 
steady  state.  The  chemistry  code  that  was  used  in  this  work 
is  based  on  the  Runge-Kutta-Mer son  algorithm  for  the  integra- 
tion of  stiff  differential  equations  originally  developed  by 
(21  ) 

Keneshea  in  his  studies  of  ionispheric  chemistry.  In  this 
scheme,  initial  CO,  02  and  He  densities  and  an  initial  elec- 
tron and  ion  density  close  to  their  expected  steady  state  val- 
ues are  specified.  A value  for  E/N  and  an  initial  vibrational 
temperature  are  also  specified  as  inputs  to  the  calculation  of 
the  electron  energy  distribution.  As  time  progresses  and  the 
CO  vibrational  temperature  changes,  the  electron  distribution 
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function  is  recomputed  and  new  rate  coefficients  calculated 
to  reflect  these  changes. 


I A list  of  the  seventy-nine  chemical  reactions  included 

in  this  calculation  is  presented  in  Table  2.  Only  a rela- 
tive few,  however,  are  really  important  and  will  be  discussed 
here.  The  rate  coefficients  for  the  reactions  and  references 
to  them  are  also  given  in  Table  2.  The  rate  coefficients 
for  the  first  twenty  reactions,  except  number  (8) , are  com- 
puted using  the  electron  energy  distribution  code  and  the 
cross  sections  discussed  above.  The  rate  coefficient  used 
for  ionization  of  vibrationally  excited  CO  (reaction  8)  is 
taken  to  be  the  same  as  that  for  ionization  of  ground  state 
CO  (reaction  5) . The  rate  coefficients  shown  in  Table  2 
for  reactions  (1)  through  (20)  were  computed  using  an  E/N  = 

1 x 10  ^ V-cm^  and  a CO  vibrational  temperature  of  3000°K 
for  a gas  mixture  corresponding  to  1.4  torr  CO,  50  m-torr  O2 
and  18  torr  He.  These  values  are  shown  for  illustrative  pur- 
poses as  they  are  typical  of  the  conditions  being  discussed 

_ 

in  this  work.  The  electron  reduced  mean  energy  (ur)  under  these 
conditions  is  about  1.0  eV.  The  gas  temperature  has  been 
taken  to  be  300°K.  The  rate  coefficients  in  reactions  (21) 
through  (79) , where  they  are  dependent  upon  electron  or  gas 
temperature,  have  been  computed  using  these  values. 

To  illustrate  how  the  parameters  critical  to  the  calcula- 
tions being  presented  here  vary  with  E/N  and  vibrational  tem- 
peratures, we  have,  in  Figures  2 through  4 plots  of  the 
ionization  rate  coefficient  of  CO,  mean  electron  energy,  and 
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TABLE  2.  Reactions  Included  in  the  Plasma 
Chemistry  Calculations. 


No.  Reactions  k Ref. 


1. 

CO 

+ 

e 

= 

cot 

+ 

e 

9.4 

(-9) 

2. 

cot 

+ 

e 

= 

CO 

+ 

e 

2.2 

(-8) 

3. 

CO 

+ 

e 

= 

CO* 

+ 

e 

6.8 

(-11) 

4. 

CO* 

+ 

e 

= 

CO 

+ 

e 

1.6 

(-8) 

5 _ 

CO 

+ 

e 

= 

co+ 

+ 

e 

+ 

e 

2.4 

(-14) 

6. 

CO 

+ 

e 

= 

c 

+ 

0~ 

2.8 

(-14) 

7. 

CO 

+ 

e 

= 

c+ 

+ 

O' 

+ 

e 

9.3 

(-17) 

8. 

cot 

+ 

e 

= 

co+ 

+ 

e 

+ 

e 

2.4 

(-14) 

9. 

°2 

+ 

e 

= 

o2t 

+ 

e 

2.3 

(-10) 

10. 

o2t 

+ 

e 

= 

°2 

+ 

e 

4.2 

(-10) 

11. 

°2 

+ 

e 

= 

°2* 

+ 

e 

] .9 

(-10) 

12. 

0 * 
U2 

+ 

e 

= 

°2 

+ 

e 

5.4 

(-10) 

13. 

°2 

+ 

e 

= 

V 

+ 

e 

+ 

e 

6 .6 

(-14) 

14. 

°2 

+ 

e 

= 

0 

+ 

0 

+ 

e 

5.1 

(-11) 

15. 

°2 

+ 

e 

= 

0 

+ 

o' 

5.5 

(-12) 

16. 

o2 

+ 

e 

= 

0+ 

+ 

o' 

+ 

e 

6.1 

(-18) 

17. 

°2 

+ 

e 

= 

0 

+ 

0* 

+ 

e 

2.6 

(-11) 

18. 

0 

+ 

e 

= 

0* 

+ 

e 

6.1 

(-10) 

19. 

0 

+ 

e 

= 

0+ 

+ 

e 

+ 

e 

4.9 

(-14) 

20. 

02* 

+ 

e 

= 

°2+ 

+ 

e 

+ 

e 

1.6 

(-13) 

21. 

0~ 

+ 

CO 

= 

CO  2 

+ 

e 

7.3 

(-10) 

22 

22. 

CO* 

+ 

CO 

= 

cot 

+ 

cot 

9.9 

(-11) 

23 

23. 

CO* 

+ 

CO 

= 

C 

+ 

CO  2 

1.1 

(-11) 

23 

24. 

co+ 

+ 

CO 

+ CO  = 

C2°2+ 

+ 

CO 

1.4 

(-28) 

24-27 

25. 

C2°2 

+ 

CO 

= 

co+ 

+ 

CO 

+ 

CO 

2.1 

(-12) 

24-27 

26. 

C2°2 

+ 

e 

= 

CO 

+ 

CO 

7.4 

(-8) 

28 

27. 

°2+ 

+ 

e 

= 

0 

+ 

0 

1.0 

(-8) 

29 

ro 

00 

°2+ 

+ 

e 

= 

0 

+ 

0* 

1.0 

(-8) 

29 

29. 

0* 

+ 

°2 

= 

0 

+ 

°2 

7.5 

(-11) 

30 

30. 

°2 

+ 

0 

+ M - 

°3 

+ 

M 

6.4 

(-34) 

31 

29 


TABLE  2.  (continued) 


No . 

Reactions 

k 

Ref . 

31. 

0* 

+ 

°3 

= 

°2 

+ 

°2 

3.8  (-12) 

30 

32. 

0 

+ 

°3 

= 

°2 

+ 

°2 

2.0  (-14) 

31 

33. 

°2 

+ 

°3 

= 

°2 

+ 

o0  +0 

1.5  (-13) 

31 

34. 

0* 

+ 

°2 

= 

4* 

C\J 

o 

+ 

0 

4.0  (-11) 

32 

35. 

+ 

04  1 

O 

+ 

e 

4" 

4-  M = 

°2 

+ 

M 

8.0  (-29) 

21 

36. 

0 

+ 

°2 

— 

0 

+ 

°2 

1.0  (-7) 

33 

37. 

0" 

+ 

0 

= 

°2 

+ 

e 

2.0  (-10) 

33 

38. 

0 

+ 

°2 

= 

°3 

+ e 

5.0 

(-15) 

33 

39. 

°3 

+ 

e 

= 

0 

+ 02 

1.0 

(-11) 

33 

40. 

0" 

+ 

°2* 

= 

°3 

+ e 

3.0 

(-10) 

33 

_ 

■4* 

41. 

0 

+ 

°2 

+ 

M 

= 

°3 

+ M 

1.9 

(-27) 

est . 

42. 

°2 

+ 

e 

+ 

M 

= 

°2~ 

+ M 

1 .0 

(-33) 

34 

43. 

0 

+ 

e 

+ 

M 

= 

0" 

+ M 

1.0 

(-31) 

21 

44. 

C 

+ 

°2 

= 

cot 

+ 0* 

3.3 

(-11) 

35 

45. 

C 

+ 

CO 

= 

CO 

+ CO 

7.0 

(-19) 

35 

46. 

C 

+ 

4 

= 

co+ 

+ 0 

2.3 

(-10) 

35 

47. 

C+ 

+ 

C02 

= 

co+ 

+ CO 

1.8 

(-9) 

32 

48. 

C+ 

+ 

o2 

= 

co+ 

+ 0 

1.1 

(-9) 

32 

49. 

co+ 

+ 

e 

= 

c 

+ 0 

9.2 

(-8) 

35 

50. 

co+ 

+ 

e 

+ 

M 

= 

CO 

+ M 

8.5 

(-27) 

33 

51. 

c 

+ 

0 

+ 

M 

= 

CO 

+ M 

1.0 

(-32) 

35 

52. 

C0+ 

+ 

°2 

= 

CO 

+ 0 + 

2.0 

(-10) 

32 

53. 

co+ 

+ 

0 

= 

CO 

+ 0 

1.4 

(-10) 

33 

54. 

C 02+ 

+ 

°2 

= 

CO 

+ CO  + o2+ 

2.0 

(-10) 

est . 

55. 

0* 

+ 

CO 

= 

°2+ 

+ CO 

1.1 

(-9) 

33 

56. 

CO* 

+ 

°2 

= 

CO 

+ 02* 

2.0 

(-10) 

23 

57. 

0* 

+ 

CO 

= 

0 

+ cot 

1.7 

(-11) 

36 

58. 

0* 

+ 

CO 

= 

0 

+ CO 

5.6 

(-11) 

36 

59. 

0 

+ 

cot 

= 

0 

+ CO 

5.9 

(-15) 

37 

60. 

°2" 

+ 

0 

= 

°3 

+ e 

5.0 

(-10) 

32 

I 
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TABLE  2.  (continued) 


No.  Reactions  k Ref. 


61. 

°2 

+ 

0 

= 

°2 

+ 

0 

3.3 

(-10) 

33 

62. 

V 

+ 

°2+ 

= 

°2 

+ 

°2 

4.2 

(-7) 

33 

63. 

°2~ 

+ 

°2 

= 

°2 

+ 

°2 

+ 

e 

2.0 

(-18) 

33 

64. 

°2 

+ 

°3 

= 

°2 

+ 

°3_ 

4.0 

(-10) 

33 

65. 

°2~ 

+ 

°2* 

= 

°2 

+ 

°2 

+ 

e 

2.0 

(-10) 

33 

66. 

0 

+ 

°3 

= 

0 

+ 

°3~ 

1.0 

(-9) 

33 

67. 

°2+ 

+ 

°3_ 

= 

°3 

+ 

0 

+ 

0 

1.0 

(-7) 

38 

68. 

°2+ 

+ 

V 

= 

°3 

+ 

°2 

2 .0 

(-7) 

38 

69. 

°3~ 

+ 

0 

= 

°2~ 

+ 

°2 

1.0 

(-10) 

38 

70. 

°3~ 

+ 

0 

= 

°2 

+ 

°2 

+ 

e 

1.0 

(-13) 

38 

71. 

°2+ 

+ 

°2 

+ 

°2 

= 

°4+ 

+ 

°2 

2.8 

(-30) 

33 

72. 

o 

ro 

+ 

+ 

°2 

+ 

M 

= 

°4+ 

+ 

M 

1.0 

(-31) 

33 

73. 

v 

+ 

°2 

= 

°2+ 

+ 

°2 

+ 

°2 

2.0 

(-13) 

33 

74. 

C 

+ 

0 

= 

°2+ 

+ 

°3 

3.0 

(-10) 

33 

75. 

cot 

+ 

He 

= 

CO 

+ 

He 

1.7 

(-17) 

39 

76. 

cot 

= 

CO 

+ 

hv 

3.4 

(+1) 

40 

77. 

°2+ 

+ 

He 

= 

°2 

+ 

He 

1.6 

(-15) 

41 

78. 

»4+ 

+ 

e 

= 

°2 

+ 

°2 

1.1 

(-7) 

29 

79. 

0+ 

+ 

e 

+ 

M 

= 

0 

+ 

M 

2.0 

(-27) 

33 

-12 

Note  that  3.8  (-12)  = 3.8  X !10  cc/molecule-sec . Also  t refers  to 
vibrational  excitation  while  * refers  to  electronic  excitation. 
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FIGURE  3:  Mean  electron  energy  as  a function  of  E/N  and  CO 

vibrational  temperature  for  a mixture  of  7%  CO/93%  He. 
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FIGURE  A:  Recombination  rate  coefficients  as  a function  o 

E/N  for  two  CO  vibrational  temperatures  (solid  line  i 
300°K,  dotted  line  is  4000°K). 


recombination  rate  coefficients. 


These 


C0+,  C202+  and  C>2+ 
were  computed  for  mixtures  of  1.4  torr  CO,  50  m-torr  02  and 
18  torr  He  corresponding  to  the  experimental  conditions  of 
Keren,  et  al.  (1°  > -*-9)  The  rate  coefficient  for  recombination  of 
C2O2  was  obtained  from  center'  J who  assumed  that  in 
his  high  pressure  experiments  the  dominant  ion  was  C202+ 
or  a higher  order  polymer  ion. 

The  results  of  the  calculation  of  the  time  evolution  of 
chemical  species  in  a CO-C>2-He  discharge  are  displayed  in 

Figures  5 through  8 . in  these  calculations  , the  electron  den- 

9-3 

sity  was  taken  to  be  5 x 10  cm  initially,  with  a gas  mixture 
of  1.4  torr  CO,  5 0 p 02  and  18  torr  He  and  E/N  equal  to  1.05  x 
10  ^ V-cm2.  The  initial  CO  vibrational  temperature  was  chosen 
to  be  1000°K  and  as  it  increased  the  electron  impact  rate  co- 
efficients were  recomputed  at  T = 1500°,  2000°,  2400°,  etc. 
The  species  CO,  02,  O,  O^*  and  He  were  included  in  the  calcu- 
lation of  the  electron  energy  distribution  function,  but  the 
presence  of  oxygen  in  small  concentrations  has  little  effect 
upon  the  shape  of  the  distribution  function,  even  in  the  high 
energy  tail. 

The  above  calculations  show  that  02+  becomes  the  dominant 
ion  in  CO  lasers  containing  02  due  to  charge  transfer  reac- 
tions and  that,  as  a consequence  of  recombination 
kinetics,  it  is  responsible  for  the  lowering  of  the  E/N 
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FIGURE  6 : CO/C^/Ee  kinetics  (cont'd). 
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FIGURE  8 : CO/O « /He  kinetics  (cont'd).  Note  that 

- - + - Z 7 

02>  0 , C , Oo  have  densities  below  1 /cc. 
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Ionization  from  the  metastable,  which  appears  in  large 
concentrations  due  to  reactions  (11)  and  (56) , 

0 2 + e C>2  + o (11) 

CO*  + 02  -*•  CO  + O/  , (56) 

has  been  included  in  these  calculations  but  is  unimportant  as 
a source  of  0+  in  comparison  to  the  charge  transfer  process. 

The  relative  rates  of  O*  ionization  and  CO  ionization  for  a 
typical  CW  flow  laser  contact  time  of  0.1  sec  are  shown  in 
Table  3.  By  reference  to  Table  3 we  see  that  the  direct 

ionization  of  CO  [reaction  (5)]  is  over  an  order  of  magnitude 

* 

larger  than  the  lonizatior  ^ate  of  metastable  0^  [reaction 

(20)].  This  observation  is  in  disagreement  with  Keren,  et 
(18  191 

al v ’ 'who  postulated  that  ionization  via  metastable  C>2  is 
the  dominant  process  for  0+  production. 

This  replacement  of  CO+  and  C202+  by  02+  as  the  major 
ionic  species  has  an  important  effect  upon  the  voltage-current 
characteristics  of  the  discharge.  For  simplicity,  if  one  assumes  that 
the  only  formation  and  removal  processes  for  electrons  are  ioni  - 
zation  and  homogeneous  recombination,  respectively,  the  elec- 
trons are  in  a steady  state  density  which  can  be  computed 
from  (recognizing  that  02+  , CO+  and  C202+  recombine  through  a 
two-body  process  at  the  pressures  of  interest  here) , 
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in  a current  regulated  discharge.  This  is  in  direct  opposi- 

/ 1 Q 1 Q\ 

tion  to  the  explanation  advanced  by  Keren,  et  al.  ’ ’ This 

lowering  of  E/N  allows  for  more  efficient  pumping  of  the  CO 
vibrational  levels  as  well  as  reduced  plasma  heating  as  a re- 
sult of  electronic  excitation  of  CO  followed  by  rapid  quench- 
ing. Both  factors  lead  to  increased  power  output.  The  cal- 
calations  also  demonstrate  that  0 atoms,  formed  predominantly 
by  electron  impact  dissociation  of  O^ , have  an  adverse  effect 
upon  laser  performance  when  O^  is  added  in  large  enough  con- 
centrations. This  is  due  to  relaxation  of  the  CO  laser  levels 
by  vibrational  ti anslational  energy  exchange  in  CO-O  collisions. 
The  details  of  these  findings,  and  others,  are  presented  in  the 
following  sections. 


IV.  RESULTS . 

A.  Ion  Kinetics  of  the  CO-C^-He  Discharge. 

It  was  pointed  out  above  that  0+  is  found  to  be  the 

dominant  positive  ion  in  CO  discharges  containing  oxygen,  even 
when  the  00  is  present  in  partial  fractions  less  than  1%.  The 
results  of  the  calculation  (Figure  7 ) are  in  agreement  with 
this  observation.  The  dominance  of  0 ^ is  due  to  O2  having 
the  lowest  ionization  potential  of  the  major  species  in  the 
discharge  and  is  formed  almost  entirely  through  charge  trans- 
fer reactions.  ^~2^2  ' is  the  dominant  ion  in  CO  dis- 

charges without  oxygen,  and  CO  botn  readily  charge  transfer 


TABLE  3 


Rate  constants  and  rates  of  individual  reactions 
at  t = 0.1  sec  using  the  reactions  in  Table  2. 


k 

Rate 

k 

Rate 

1. 

9 . 2 2 ' 9 , 

5. 36 18 

26. 

7 . 3 2 " ° 8 , 

2. 

2 . 2 4 " 8 , 

5 . 1 91 8 

27. 

1.01'08, 

1 . 30 13 

3. 

~ . 67  " 1 1 , 

t /ir16 

4.46 

28. 

i.or08, 

1.3013 

4 . 

1 . 57  " 8 , 

1 . 4 ? 1 3 

29. 

7.50'11  , 

1 . 6 4 1 4 

5. 

2.98"14, 

1 *7  ’ ^ ^ 

x . / O 

30. 

6. 45" 34  , 

i - n 1 5 

1 . o 7 

6. 

5.29-14, 

1 

1.9i 

31  . 

- 1 9 

3.3  0 , 

3 . 0 1 1 0 

7. 

1 . 1 2 " 1 6 , 

6. 50 10 

32. 

- 1 4 

2.00  , 

C.1412 

8. 

2.9S*14  , 

6 . 9 0 1 

33. 

1 . 50" 1 

9. 79  12 

9. 

2.3S"10, 

1 . 1 6 1 3 

34  . 

4.00"11 , 

1 ] 

2.64 

10. 

4.16'10, 

3.941-5 

35. 

8 . 00*  29  , 

T ,-10 

3 . 1 o 

11  . 

1 .88“ 10 , 

9 . 3 5 1 4 

36. 

1.00'C/ , 

■>  n 

1 .471" 

12. 

5.44'10, 

2.9515 

37. 

2 . 0 0 " 1 0 , 

c --11 
5 . .5.9 

13. 

8 . 0 2 " 1 4 , 

3.9S11 

58. 

5. 00* 15  , 

2.S400 

14  . 

5. Si'11  , 

2 . 7 3 1 4 

39  . 

1 .00" 11 , 

1 1 

1.79* 

15. 

5 . 9 7 " 1 2 , 

2 . 9 31 ^ 

40. 

3.00"10, 

1.3b11 

16. 

7 . 4 6 " 1 8 , 

3 . 7 0 7 

41  . 

-27 

1.90  , 

8. 5 3° 7 

17. 

2.99']1  , 

1 . 4 8 1 4 

42. 

1.00'33, 

1 .52°  9 

18. 

6 . 16" 10  , 

i ■ ,16 

1 . “4  4 

43. 

1 . oo'^1 , 

7 . 1 5 1 1 

19. 

5 . 9 7 " ■*  4 , 

1 7 

1.39 

44  . 

s . J>0  , 

4.3815 

20. 

1 .90* 13 , 

1 9 

1 .03 

45. 

- 1 9 

6.98  , 

2 . 2 3°  8 

21. 

7.30"1C, 

4 . 8 7 1 3 

46. 

2.30  1 ° , 

7.8712 

22. 

9 . 90*  8 , 

3 . 9 5 1 6 

47  . 

1.S0"09, 

5.4210 

23. 

•s  , n - 11 

i . i0  , 

4 . 3 9 1 5 

48. 

1 . 1 0 " 0 9 , 

1 . 03 1 0 

24. 

- V 

1.43  " , 

1 2 

8 . 4 : 1 L 

49. 

9 . 2 0 " ° 8 , 

8.05 

25. 

.19 

2.10  , 

4 . 44  1 4 

50. 

S . 4 0 ' fc  ''  , 

2 . 2 5 1 9 
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TABLE  3 (cont'd) 


51 . 

- 3? 

1.00 

1 . 91 1 “ 

52. 

2.00'10, 

1 .08 10 

53. 

1 .40' 10, 

2.  1013 

54. 

: . oo' iU , 

, .,12 
O . 0 1 

55. 

1 . 10' 09 , 

7 1^-'^ 

c . - i 

56. 

2 . 0 0 ' 1 0 , 

6 . 8 2 1 4 

57. 

1 . 72'  , 

4.4215 

58. 

5 . 58  * 1 1 , 

1.4316 

5S . 

5 . 9 0 ' 1 5 , 

2.3416 

60. 

5 . 0 0 " 1 0 , 

i.S309 

61 . 

3 . 2 8 ' 1 0 , 

1 . 2009 

62. 

4 . 2C'07, 

8 . 4207 

63. 

2.00'18, 

1.56°° 

64  . 

4.00'iU, 

1 . 1 2 0 6 

6 5 . 

2.00" l0, 

1 . 69  08 

66  . 

i.oo'09, 

2 . 0 5 ° 9 

67. 

1 . oo' 0/ , 

9 . 6 7 0 7 

OO 

o 

2 . 0 0 ' ° ' , 

0 8 

1.93 

60. 

i.oo'10. 

, -.09 
1 . / 0 

70. 

1 . oo'i3, 

1.70?6 

7 1 . 

2 . 8 0*  30  , 

1.7 5° 9 

— 7 

1 -0«  , 

1 . 1 3 1 1 

7 3 . 

2 . 00' 1 

7 . OO07 

7 4 . 

- -,.-10 
j . 0 0 , 

1 . 4 1 1 1 

7 5 . 

1 . 

. 70 

-17 

> 

3 . 

O 

o 

76. 

3, 

. 36 

+ 01 

t 

2. 

.ll17 

77. 

1 . 

.60 

-15 

y 

1. 

, 1 9 1 3 

78. 

1 . 

.11 

-07 

y 

3. 

. os09 

79. 

2 . 

.00 

-27 

y 

1. 

ll09 

k . [e 
ion 


- k , [e 

recomb 


If  the  ion  and  electron  densities  are  assumed  to  be  equal 
[i.e.,  one  ion  dominates  and  negative  ions  are  unimportant 
at  steady  state . 


= 0 ->  [e  ] 


k . 
ion 

"recomb 


As  can  be  seen  in  Figure  4,  the  rate  coefficient  for  recom- 
bination of  02+  is  substantially  smaller  than  that  of  either 
C0+  or  C202+  • This  results  in  an  increased  electron  density 
when  02+  dominates  so  that,  under  constant  current  operation, 
the  electric  field  must  be  reduced.  A lower  electric  field  in- 
creases the  efficiency  of  electron  impact  excitation  of  the 
CO  vibrational  states  as  well  as  reducing  the  electron  energy 
channeling  into  electronic  excitation  and  hence  heat,  event- 
ually resulting  in  increased  power  output. 


As  a test  of  this  hypothesis  , a simple  model  calculation  was 
performed  in  which,  given  the  recombination  rate  coefficient 
(as  a function  of  electron  temperature) , E/N  and  CO  vibra- 
tional temperatures  were  varied  and  the  ionization  rate  co- 
efficient calculated  in  order  to  achieve  a chosen  current  den- 

2 

sity.  Current  density,  j (amp/cm  ) is  related  to  electron 

-3 

drift  velocity,  v^  (cm/sec)  and  number  density,  Nc  (cm  ) by 
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where  e = 1.6x10  coul.  and  where  the  drift  velocity 
can  be  computed  from  the  electron  energy  distribution  func 
tion. 


If  a similar  steady  state  analysis  is  carried  out  for 
excitation  of  CO(v=l)  by  electrons  and  de-excitation  by  super- 
* elastic  collisions  with  electrons,  which  is  predicted  (see 

Table  3)  to  be  one  of  the  major  loss  processes  for  low  vi- 
brational levels,  we  find  at  steady  state 

[CO ( v=l ) 1 _ k0+l 

I [CO  (v=0)  ] ' 

[ Using  these  two  expressions  and  the  equation  for  current  den- 

sity we  can  choose  an  E/N  and  CO  vibrational  temperature,  solve 
the  Boltzmann  equation  obtaining  rate  coefficients,  temperature 
and  drift  velocity,  and  --alculate  a current  density  and  new  vi- 
brational temperature.  is  process  can  be  repeated  until  the 
assumed  and  computed  vibrational  temperatures  are  equal  and 
. the  desired  current  density  has  been  obtained.  This  approach 

is  simple,  but  contains  the  essence  of  the  processes  involved 
in  establishing  the  electron  density  in  the  plasma.  The  great- 
est inaccuracy  lies  in  the  assumption  concerning  vibrational 
excitation  and  de-excitation  of  CO(v=l)  since  there  are  other 
important  processes,  such  as  excitation  and  de-excitation  to 
and  from  CO(v=2,3,  . . . ),  that  affect  the  population  of  the 
v=l  level. 

I 

I 
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whore  e = 1.6  x 10  coul.  and  where  the  drift  velocity 
can  be  computed  from  the  electron  energy  distribution  func- 
tion. 


If  a similar  steady  state  analysis  is  carried  out  for 
excitation  of  CO(v=l)  by  electrons  and  de-excitation  by  super- 
elastic collisions  with  electrons,  which  is  predicted  (see 
Table  3)  to  be  one  of  the  major  loss  processes  for  low  vi- 
brational levels,  we  find  at  steady  state 

[CO  (v=l)  ] = k0-»l 

[co (v=0)  ] ' klH>0  • 

Using  these  two  expressions  and  the  equation  for  current  den- 
sity we  can  choose  an  E/N  and  CO  vibrational  temperature,  solve 
the  Boltzmann  equation  obtaining  rate  coefficients,  temperature 
and  drift  velocity,  and  calculate  a current  density  and  new  vi- 
brational temperature.  This  process  can  be  repeated  until  the 
assumed  and  computed  vibrational  temperatures  are  equal  and 
the  desired  current  density  has  been  obtained.  This  approach 
is  simple,  but  contains  the  essence  of  the  processes  involved 
in  establishing  the  electron  density  in  the  plasma.  The  great- 
est inaccuracy  lies  in  the  assumption  concerning  vibrational 
excitation  and  de-excitation  of  CO(v=l)  since  there  are  other 
important  processes,  such  as  excitation  and  de-excitation  to 
and  from  CO(v=2,3,  . . . ),  that  affect  the  population  of  the 
v=l  level. 
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Despite  its  defects,  this  calculation  gives  results  rea- 


1 


sonable  enough  to  illustrate  the  point  concerning  the  rela- 
tionship between  recombination  rate,  E/N  and  efficiency  of 
vibrational  excitation.  The  results  are  presented  in  Table 
4.  To  achieve  similar  current  densities  and  vibrational 
temperatures  using  the  0 ^ recombination  coefficient  required 
a 10%  decrease  in  E/N  from  the  value  needed  using  the  recom- 
bination coefficient  for  CO+ . The  lowering  of  E/N  resulted 
in  a 34%  decrease  in  the  fraction  of  energy  flowing  into  CO 
electronic  states  and  a 9%  increase  in  the  fraction  being 
channeled  into  v=l  alone.  Thus,  in  contrast  to  Keren,  et.  al , 
we  have  shown  that  O + directly  affects  both  the  E/N  of  the 
plasma  and  the  excitation  of  CO  vibrational  energy,  in  agree- 
ment with  the  experimental  observations. 

B.  Negative  Ions  in  a CO  Discharge. 

The  time  development  of  the  negative  ion  densities  is 
shown  in  Figure  8.  Although  O is  rapidly  formed  by  disso- 
ciative attachment  of  CO  and  reactions  (6)  and  (15) , it 

is  rapidly  removed  by  detachment  collisions  with  CO  [reaction 

(21)].  This  is  in  contrast  to  the  chemistry  of  CO2  lasers  as 

(42) 

discussed  by  Nighan  and  Wiegand  'and  by  Garscadden  and  his 

(A3  44)  , 

collaborators,  ’ where  negative  10ns  appear  to  control  the 
stability  of  the  discharge  at  low  E/N.  The  results  obtained 
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TABLE  4. 

Results  of 
ionization- 

the  simplified 
-recombination 

! 

calculation  of 
processes . 

U s i i 

Us  inn  kn’- 

F/M  (V-cn“) 

1 . 1 2 X l 0" ] 5 

i .oixi y 10 

[c  ] ( cn  '"’j 

4 . 93x1  0 + 1(1 

4 . 99x1 0+1° 

i (amp/ cm*) 

2 3 . 4 x 1 0 ‘ 

n ivi n"J 

e (eV) 

] . 27 

l . oy 

vdriFt  (C’:1/SC^ 

-i  - i n + G 

2 . 9 / x i 0 

+ 6 

2.78x1': 

rrco-) 

2i.3  : 

14.9  0 

r<cfW 

1 S . 1 7 

19.8  S, 

here  are  consistent  with  the  observations 


that  additions  of 


(42) 

CO  to  CO2  lasers  stabilize  the  discharge  through  electron  de- 
tachment reactions  with  the  negative  ionic  species. 

In  order  for  attachment  to  Fe(CQ)5  impurity  to  be 
important  as  pointed  out  by  Center^2®)  the  impurity  level 
would  have  to  be  greater  than  145  ppm.  This- level  is 
substantially  above  the  maximum  level  found  by  Center  in 
his  experiments,  therefore,  we  expect  that  impurity  attach- 
ment is  not  important  to  the  conclusions  reached  in  this  study. 

C.  The  Role  of  Atomic  Oxygen  in  the  Discharge. 

As  mentioned  earlier  in  this  section,  at  higher  par- 
tial pressures  the  laser  output  power  is  found  to  decrease  and 
the  neutral  temperature  increase.  Our  calculations  indicate 
that  both  of  these  effects  are  due  to  the  formation  of  0 atoms 
in  the  plasma.  The  cross  sections  for  electron  impact  excita- 
tion of  predissociative  O2  electronic  states  are  large.  This 
leads  to  a large  production  of  oxygen  atoms  (reactions  14  and 
17)  and  a build-up  in  the  0 atom  concentration  in  the  dis- 
charge. The  calculated  time  development  of  the  several  neu- 
tral oxygen  species  is  shown  in  Figure  6. 

For  very  small  concentrations  of  O2  as  an  additive  to  CO  la- 
sers the  dissociation  into  0 atoms  has  little  effect  upon  la- 
ser performance.  Above  some  threshold  0 ^ concentration,  the 
0 atom  density  increases  to  a level  where  vibrational-trans- 
lational  energy  transfer  collisions  between  CO  and  0 (reac- 
tion 59)  is  predicted  to  become  the  dominant  de-excitation 
mechanism  among  the  CO  vibrational  levels  responsible  for  la- 
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ser  action.  Since  it  is  the  hiqhcr  vibrational  levels  that 
dominate  the  laser  output  spectrum,  the  presence  of  0 atoms 
will  degrade  the  laser  .output . A sample  of  available  vibra- 
tional de-excitation  rate  coefficients  for  CO(v=l)  and  C0(v=12) 
are  shown  in  Table  5.  For  the  low  vibrational  states  (v=l 
to  v=8)  superelastic  collisions  with  electrons  is  the  primary 
de-excitation  process.  In  the  absence  of  oxygen,  radiation 
(reaction  76)  and  VT  collisions  with  He  (reaction  75)  are  pri- 
marily responsible  for  CO  vibrational  de-excitation.  The  ra- 
diative rate  for  C0(v=12)  molecule  is  given  by  the  Einstein  A 
coefficient  in  Table  4.  The  deactivation  rate  for  C0(v=12) 
by  collision  with  He  is  ky^tHe]  = 262  sec  ^ for  18  torr  of  He. 
This  is  approximately  the  same  as  the  radiation  rate.  Using 
a superelastic  rate  coefficient  for  the  v = 12-*-ll  transition 
equal  to  that  for  the  v = l-*0  and  an  electron  density  of  10'*'^ 
gives  a de-excitation  rate  of  10  -100  sec  Hence  radiation 

and  He  V-T  exchange  are  likely  to  be  the  dominant  de-excitation 
mechanisms. 


However,  when  O atoms  are  present,  even  in  small  concen- 
trations, this  may  strongly  affect  vibrational  relaxation 
since  the  CO-O  V-T  rate  coefficient  for  the  lowest  vibrational 
levels  is  more  than  one  hundred  times  larger  than  that  for 
CO-He/^^f  the  same  scaling  with  v is  assumed  for  V-T  relaxa- 
tion by  O atoms  as  exists  for  the  He  V-T  process  (a  substan- 
tial underestimate  based  on  SSH  mass  effect) , an  O atom  par- 
tial pressure  of  only  52  m-torr  would  be  needed  for  CO-O  V-T 
to  be  equal  to  the  radiation  or  the  CO-He  V-T  process.  Laser 
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TABLE  5 


De-excitation  of  CO  vibrational  levels. 


r 

Process  v 1 ->  0 v = 1 2-»  1 1 

;1t  40  - ] - i 

Radiation  A-  33.6  sec  239.9  see 

| ^ 

39  -t  — 

V-T  with  He  , k=  I.7xl0']7  cm^/scc  9 . Ox  1 0 " 1 6 cm'Vsec 

fT  = 300°K 1 
t ' gas 

.1 

37  i-  ~J 

V-T  with  0 , " 5.9xlO‘1S  " -2.9X10"1*' 


Svipere  las  tic 
with  c , 


'live  value  shown  i s computed  as  sun  in.  the  same  scaling 
with  v as  the  CO-ile  V-T  process. 

X X 

The  rate  coefficient  is  typical  iy  i n this  range  wi  tli  its 
exact  value  depending  upon  the  electron  energy  distribution 

f unc t i on . 


performance  will  be  degraded  when  the  contribution  of  the 
CO-O  V-T  rates  is  non-negligibie  in  comparison  to  the  fixed 
radiation  and  CO-IIe  V-T  rates.  The  scaling  of  vibrational 
rate  coefficients  with  increasing  vibrational  level  is,  in 
general,  uncertain.  However,  the  collision  system  of  0(^P)  + 
CO(1Z,v'),  which  correlates  with  the  triplet  state  of  CO 2, 
offers  the  possibility  of  a strong  "chemical"  interaction 
through  curve  crossing  to  the  state  and  enhanced  en- 

ergy transfer  possibilities  associated  with  a short  lived  tri- 
atomic  complex. 

Although  the  above  arguments  are  somewhat  qualitative  due 
to  a lack  of  reliable  rate  constant  data  information,  our  chem- 
istry calculations  predict  substantial  dissociation,  larger 
than  20%,  for  time  scales  of  the  order  of  flow  contact  times 
in  typical  CO  lasers.  For  dissociation  levels  of  this  order,  0 
atom  VT  processes  will  substantially  reduce  the  vibrational 
content  of  CO,  leading  to  reduced  laser  output  as  well  as  in- 
creased plasma  heating. 

The  effect  of  higher  C>2  additions  on  the  shape  of  the 

experimentally  observed  CO  vibrational  distribution  up  to 

v = 30  also  clearly  shows  the  emergence  of  a strong  VT 

(45  46) 

process,  attributable  to  0 atoms.  ’ 


l 
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It  can  be  seen  from  the  rate  coefficients  in  Table  5 


that  vibrational  relaxation  by  0 2 will  be  unimportant  for  the 
concentrations  being  discussed  here.  The  same  will  be  true 
for  small  concentrations  of  CO^ , which  is  formed  primarily  by 
reactions  (21)  and  (23)  : 


+ CO  ->• 

C02  + e 

(21) 

+ CO  -V 

c + co2 

(23) 

The  density  of  C02  is  still  increasing,  however  (Figure  7) , 
at  the  end  of  the  calculation  indicating  that  there  may  be  a 
large  buildup  of  this  species  in  closed  CO  systems.  In  flow- 
ing discharges  with  a residence  time  of  .1  - .5  sec  this  may 
not  be  a problem,  but  in  a sealed  system,  C02  may  be  present 
in  large  enough  concentrations  to  be  a major  source  of  vibra- 
tional de-excitation. 


D.  Additional  Discussion  of  CO  Laser  Chemistry. 


We  propose  that  carbon  in  the  CO  laser  is  formed  pri- 
marily as  a product  in  the  quenching  of  electronically  ex- 
cited CO,  i.e.,  reaction  (23).  Reactions  (22)  and  (23), 


★ 

CO 

+ CO  ->• 

•f*  + 

CO  + CO 

(22) 

CO* 

+ CO  -> 

c + co2  , 

(23) 

★ 

correspond  to  two  channels  for  de-excitation  of  CO  by  CO. 

(23) 

The  total  rate  coefficient  is  known;  but  the  branching  ratio 
is  not.  It  has  been  assumed  in  these  calculations  that  10% 
of  the  CO  + CO  reactions  follow  the  C + CO2  branch.  The  C 
production  rate  is  critically  dependent  upon  this  branching 
ratio  and  there  is,  unfortunately,  no  information  on  the 
branching  ratio  or  even  on  the  carbon  concentrations  found  in 
CO  discharges.  The  carbon  is  readily  removed  by  reaction  (44) 

C + 02  ->  COf  + 0*  (44) 

with  02-  This  is  a satisfactory  explanation  of  the  observa- 
tion that  the  addition  of  02  to  a CO  discharge  reduces  or 
eliminates  carbon  deposition  in  the  system.  The  products  of 
this  reaction  have  been  assumed  to  be  CO  and  0(^D)  as  these 
are  spin  allowed  and  energetically  accessible. 

Finally,  mention  should  be  made  of  the  possible  role  of 
0(^D)  in  the  discharge.  It  is  rapidly  formed  by  reactions 
(17) , (18)  and  (28) : 


5? 


T 


°2 

+ e -*■ 

0 + 0 + e 

•k  — 

(17) 

0 

+ e -*• 

0 + e 

(18) 

+ 

* 

°2 

+ e -* 

0 + 0 

(28) 

and,  perhaps,  (44)  as 

discussed  above. 

★ 

The  0 is  then  rapidly 

removed  in  reactions 

(57)  and  (58)  : 

* 

0 

+ CO  -*■ 

o 

+ 

n 

o 

a- 

(57) 

* 

0 

+ CO  -+ 

0 + CO  . 

(58) 

Although  0*  has  little  effect  upon  the  chemistry  of  the  sys- 
tem it  may,  in  this  way,  be  an  important  source  for  heating 
the  gas. 

Using  the  reaction  rates  for  both  0(^D)  quenching  into 
translational  energy  (reaction  58)  and  relaxation  of  CO1  into 
translational  energy  (reaction  59)  as  given  in  Table  3,  the 
total  energy  flux  from  these  reactions  is  sufficient  to  ex- 
plain the  observed^^  gas  temperature  increase  with  increasing 
pressure  beyond  the  optimum  pressure  for  laser  output. 

V.  HETEROGENEOUS  VERSUS  HOMOGENEOUS  ELECTRON  LOSS  PROCESSES 

In  low  pressure,  small  diameter  laser  discharges,  ambi- 
polar  diffusion  can  dominate  the  electron  loss  mechanism. 

Here  we  consider  the  competition  between  homogeneous  recombin- 
ation and  ambipolar  diffusion  in  controlling  the  electron  den- 
sity in  a CO  laser  system. 
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In  order  to  estimate  the  relative  contribution  to  the 
electron  loss  rate,  the  time  scale  associated  with  diffusion 
and  recombination  must  be  estimated.  The  time  scale  for 
diffusion  loss  is  found  from  a solution  to  the  radial  dif- 
fusion equation  as 


dif  f 


kT 


J+ 


where  y+  is  the  ion  mobility,  Tg  is  the  mean  electron 
temperature  and  A is  the  characteristic  diffusion  length 
(A  >_  R/2.4  where  R is  the  tube  radius).  For  a 20  torr, 

300 °K  plasma  in  which  CO+  is  the  dominant  ion  (taking  R = 

(47) 

1.5  cm  and  T = 0.9  eV) 
e 

CO+  c . -4 

T , • c c = 5.1  sec . 

dif  f 


The  ambipolar  diffusion  time  for  0„  as  the  dominant  ion  is 

1 

+ / ® 2 -4  \ 

only  slightly  slower  than  that  of  CO  \Tciiff  = ^'^  seel  under 

(47) 

the  same  conditions.  Evidence  suggests  that  the  diffusion  of 


ion  clusters  in  the  case  of  C2O2  is  not  substantially  differ- 
ent from  the  parent  ion. 

The  time  scale  for  homogeneous  recombination  can  be  es- 
timated from  the  rate  of  the  recombination  reactions,  (26)  in 
the  case  of  C202+  and  (27)  + (28)  in  the  case  of  02+  . Thus, 


and 


C2°2  . -4 

T recomb  = 4*6  sec 


° 2 . ?-3 

1 recomb 
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These  estimates  show  that  the  ambipolar  loss  time  scale 
and  the  homogeneous  recombination  time  scale  are  of  the  same 
order.  At  higher  pressures  and  lower  temperature  operation 

the  homogeneous  recombination  process  will  dominate  and  the 
calculations  presented  above  on  the  influence  of  oxygen  m 
the  electron  density  will  more  quantitatively  apply.  Under 
the  conditions  of  the  calculations  presented  in  this  study, 
the  quantitative  predictions  will  be  reduced,  but  the  substi- 
tution of  O*  as  the  dominant  ion  will  still  clearly  reduce 
the  rate  of  electron  loss  and  produce  a corresponding  increase 
in  electron  density  at  constant  E/N. 
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VI.  SUMMARY  AND  CONCLUSION. 

We  have  attempted  to  demonstrate  a consistant  interpre- 
tation to  the  effects  of  small  amounts  of  molecular  in  alow 
pressure  CW  CO/Ke  laser  system.  The  observation  that  very  small 
amounts  of  lead  to  a reduction  in  E/N , a reduction  in 
plasma  temperature  and  an  increase  in  output  power  we  attrib- 
ute to  the  emergence  of  0 ^ as  the  dominant  positive  ion. 

Since  0 * recombines  more  slowly  than  C0+  or  its  polymers  at 
the  E/N  of  interest,  the  same  current  can  be  maintained  by  a 
lower  E/N.  The  reduction  in  E/N  results  in  a more  favorable 
fraction  of  the  electron  energy  into  vibration  with  a corre- 
sponding decrease  of  electron  energy  into  CO  electronic  states. 
The  thermal  energy  produced  by  electronic  quenching  (E-V/T) 
is  reduced  and  the  plasma  cools.  The  simultaneous  improved 
vibrational  excitation  and  reduced  plasma  temperature  leads 
to  improved  laser  energy  output. 

At  higher  0 2 additions,  the  0 atoms  generated  primarily 
by  predissociation  of  O2  lead  to  direct  deactivation  of  the 
upper  vibrational  levels  of  CO  responsible  for  lasing  through 
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VT  processes  as  well  as  increased  plasma  temperatures  through 
CM^D)  formation  by  electron  impact  followed  by  quenching  by 
CO  through  E-V/T  processes.  These  two  0 atom  effects  lead  to 
a degradation  of  laser  output. 

The  dynamics  of  carbon  formation  and  loss  is  proposed  to 
occur  through  electronically  excited  CO  but  is  not  found  to 
be  important  in  the  laser  dynamics  except  to  influence  the 
laser  lifetime.  This  latter  effect  occurs  when  the  carbon 
deposits  increase  cavity  losses  such  that  lasing  ceases. 

The  importance  of  plasma  chemistry  processes  involving 
minor  constituents  attests  to  the  complexity  of  molecular  dis- 
charge sustained  lasers.  The  dominance  of  recombination  kin- 
etics as  opposed  to  ionization  processes  suggests  that  other 
additives  might  enhance  laser  output  in  ways  previously  un- 
appreciated . 

Lastly,  we  would  like  to  remark  on  the  role  of  super- 
elastic processes  in  the  CO  laser  system.  It  is  clear  that 
superelastic  feed-back  from  the  vibrational ly  excited  CO  is 
critical  in  determining  observed  ionization  rates,  and  hence 
electron  densities.  We  would  like  to  suggest  that  super- 
elastic processes  are  also  responsible  for  the  normal  nega- 
tive E vs.  I characteristics  measured  in  low  temperature  CO 
lasers  and  for  the  deviation  to  a positive  E vs.  I character- 
istic observed  under  higher  temperature  operation.  We  sug- 
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gest  that  without  C>2  / the  E/N  of  the  discharge  is  high  but 
that  at  liquid  N0  cooling  the  CO  vibrational  modes  are  still 
sufficiently  excited  that  the  superelastic  feed-back  deter- 
mines the  ionization  rate  and  electron  density.  As  the  E/N 
is  reduced,  a more  favorable  vibrational  excitation  by  elec- 
trons is  produced  coincident  with  a reduction  in  plasma  tem- 
perature. The  subsequent  increase  in  vibrational  excitation 
leads  to  increased  superelastic  feed-back,  greater  ionization 
and  an  increased  electron  density,  i.e.,  current.  This  be- 
havior demonstrates  a negative  E vs.  I characteristic. 

Under  very  high  plasma  temperatures  such  as  maintained 
by  Keren,  e_t  erL,  the  superelastic  contribution  is  minimal. 

Therefore,  the  ionization  is  determined  by  the  high  energy 
tail  of  the  electron  distribution  as  provided  by  the  E/N. 

When  the  E/N  is  reduced  the  high  energy  tail  is  depressed 
leading  to  reduced  ionization  rates  and  lower  electron  den- 
sities, i.e.,  a reduction  in  current.  This  behavior  is  char- 
acterized by  a positive  E vs . I characteristic. 

It  is  clear  from  this  suggestion  that  low  temperature 
operation  or  the  presence  of  additives  would  give  negative 
E vs.  I behavior  in  agreement  with  observation.  Since  the 
only  positive  E vs.  I characteristic  has  been  reported  under 
high  temperature  operation,  we  believe  this  explanation  is 
consistent  with  available  experimental  observation. 
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Abstracts  of  papers  presented  at  the  5th 
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St.  Louis,  Mo.  April  18-20,  1977 


(1) *An  Analytical  Theory  of  Vibrational  Relaxation 

for  Anharmonic  Oscillators  Under  Strongly  Pumped 
Conditions,  S.H. Lam, Princeton  University. 

(2)  Measurements  of  CO  Vibrational  Distributions  in 
Mixtures  with  N~,  0?,  A r and  He,  A.J.Lightman  and 
E.R. Fisher,  RIFS,  Wayne  State  University. 

(3)  Comparison  of  Experimental  and  Theoretical  Vibrational 
Population  Distributions  in  Electric  Discharge  CO-He 
Mixtures,  E.P. Fisher,  R I E S , Wayne  State  University, 

H. Rabitz  and  S.H. Lam,  Princeton  University. 

(4)  Plasma  Chemistry  and  Energy  Transfer  Processes  in 
CO-O^-He  and  CO-Nn-He  Mixtures,  E.P. Fisher  and  G. 
Abranam,  RIES,  Wayne  State  University. 


*Professor  Lam's  first  paper  was  supported  by  AFOSR 
(F44620-73-0059) . It  is  included  here  since  it 
forms,  together  with  the  following  3 papers,  a 
fully  consistent  interpretation  of  the  dynamics  of 
high  lying  vibrational  states  in  CO. 


An  Analytical  Theor v o f V ib r a t_i ona  1 Relaxa- 
tion for  Anharmonic  Oscillators  Under  St rongly 
Pumped  Conditions.  S.  11.  I .AM,  Princeton  U . * — 

An  analytical  theory  for  the  quasi-steady  vibra- 
tional population  distribution,  N^,  of  a dia- 
tomic gas  under  strong  vibrational  excitation 
will  be  presented,  showing  explicitly  the  de- 
pendence of  N-:  on  the  relevant  VV  and  VT  kinetic 
rates,  temperature,  anharmonicity  of  the  energy 
levels,  and  the  strength  of  the  external  pump- 
ing. It  is  shown  that  from  the  large  collections 
of  VV  and  VT  rates  needed  for  the  master  kinetic 
equations,  two  parameters  and  two  functions 
suitably  extracted  are  the  only  important  in- 
formation required.  The  location  of  the  Treanor- 
Plateau  and  the  Plateau-Boltzmann  boundaries, 
the  ground-level  vibrational  temperature,  the 
amount  of  vibrational  energy  dissipated  by  VV 
and  VT  collisions  (in  quasi-steady  situations) 
and  the  magnitude  of  the  Plateau  are  obtained. 
These  explicit  analytical  results  can  be  used  to 
deduce  quantitative  information  from  experimental 
measurements  on  the  kinetic  rates,  particular- 
ly the  VT  rates  at  high  quantum  numbers. 


*This  research  was  supported  by  AFOSR  (F44620- 
73-0059). 


Abstract  Submitted  to  the 
5*'1  Conference  on  Chemical  and  Molecular  Lasers 

April  18,  19,  20,  1977 


Measurements  of  CO  Vihrationa 1 Distributions 

in  Mixtures  with  N-?, O2 , Ar  and  He  ■ * A.  J . LICHTMAN 

and  E.R.  FISHER,  RIES,  Wayne  State  U. — The  first  and 
second  overtone  infrared  emission  from  a discharge 
sustained  CO  laser  plasma (1)  in  mixtures  with  N2,  O2 , 

Ar  and  He  has  been  measured.  A 2.5cm  diameter  fast 
flow  lase~  tube  was  liquid  nitrogen  cooled  and  was 
operated  with  CO  partial  pressures  from  0.5  to  1.5  torr 
total  pressures  f rom  4 to  16  torr  and  currents  from  5 
to  15  ma.  Least  squares  analysis  of  the  data  yielded 
the  rotational  temperature  and  the  relative  vibrational 
populations  in  CO  from  v = 2 to  v = 30.  Detailed  model 
calculations  on  the  CO-He  system  have  yielded  energy 
transfer  scaling  relationships  to  high  vibrational 
states. (2)  Under  constant  partial  pressures  and  for 
varying  currents  from  5 to  15  ma  in  CO-he , a clear 
radial  variation  in  the  electron  density  is  indicated. 
Measurements  in  C0-N2-He  mixtuies  show  the  enhancement 
in  the  distribution  due  to  N2-CO  W pumping.  Measure- 
ments on  CO-N2-AR  mixtures  show  higher  rotational 
temperatures  than  with  He  and  give  insight  into  the 
effect  of  mass  on  the  scaling  of  VT  processes.  Detailed 
measurements  on  increasing  O2  additions  to  both  CO-He 
ar.d  C0-N2~He  mixtures  show  enhanced  distributions  for 
small  Ot  additions  and  clear  VT  effects  in  the  upper 
states  for  higher  O2  additions. 

’'partially  supported  by  DARPA  through  ONR-Boston. 

(1) a.J.  Lightman  and  E.R.  Fisher,  Appl.  Phys.  Lett.  29, 
593  (1976). 

(2) e.R.  Fisher,  H.  Rabitz  and  S.H.  Lam,  following 
paper . 

0)w.L.  Morgan  and  E.R.  Fisher,  manuscript  in  pre- 
paration . 
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Comparison  oi  Experimental  and  Theoretical 
Vibrational  Population  Distribution:;  in  Elect ric 
Discharge  CO-Hc  Mixtures.*  E.  R.  FISHER,  Wayne 
S:  .a  f U.  , H.  lUiilTZ,  and  S.  H.  LAM,  Princeton  U. 
Theoretical  predictions  cf  vibrational  population 
distributions  (N^)  have  always  been  hampered  by 
the  uncertainties  of  the  relevant  kinetic  rates 
in  the  master  kinetic  equations.  The  present 
paper  compares  detailed  experimental  measurements 
of  Ni  in  CO-He  mixtures  (T°  = 140°K,  0 £ i < 30) 
with  theoretical  calculations  using  various  rate 
models.  These  calculations  are  shown  to  be  par- 
ticularly sensitive  to  the  VT  rates  assumed. 

Using  the  recent  theoretical  rates  of  Verter  and 
Rabitz,  reasonable  agreements  between  theory  and 
experiments  are  obtained.  By  slightly  adjusting 
the  Verter-Rabitz  rates,  excellent  agreements  have 
been  obtained.  The  Verter-Rabitz  rates  have  a 
distinctly  slower  quantum  number  dependence  than 
conventionally  assumed  (extrapolations  using  SSH) 
at  this  temperature,  and  the  results  of  the  present 
study  give  strong  support  to  this  distinctive 
feature.  The  effects  of  W rate  models  are  also 
studied  and  will  be  discussed. 

*This  research  was  supported  partially  by  DARPA 
through  ONR-Boston  and  partially  by  AFOSR. 


5*^'  Conference  on  Chemical  and  Molecular  Lasers 


April  18,  19,  20,  1977 


Plasma  Chemistry  and  Energy  Transfer  Processes 
in  C0-02~He  and  CO-N^-He  Mixtures.  E.R.  FISHER  and  C. 
ABRAHAM,  RIES,  Wayne  State  U. — Experimental  evidence 
as  well  as  past  model  calculations  have  demonstrated 
the  important  coupling  between  the  plasma  ion  chemistry, 
neutral  chemical  reactions  and  energy  transfer  processes 
in  affecting  the  output  and  efficiency  of  discharge 
sustained  CO  molecular  lasers.  When  a small  amount  of 
O2  is  added  to  a CO-He  mixture  the  output  power  is 
known  to  be  increased  while  larger  additions  cause 
the  output  to  degrade.  In  past  work(l)  we  have  used 
a coupled  laser  model  incorporating  a non-Boltzmann 
electron  energy  code,  a transient  ion  and  neutral 
chemistry  code,  and  a detailed  energy  transfer  model 
to  explain  these  and  other  O2  observations.  The  laser 
enhancement  is  a result  of  the  recombination  dynamics 
of  02+  compared  to  C0+  (or  C202+)  while  the  output 
degradation  is  a result  of  0-atom  VT  relaxation.  The 
plasma  chemistry  and  energy  transfer  codes  used  above 
have  been  extended  to  include  N2  and  N containing 
species.  In  particular,  relatively  high  concentrations 
of  CN,  and  the  electronically  excited  A^tt  and  b2£ 
states  of  CN,  have  been  observed  in  C0-N2~He  mixtures. 

A decrease  in  laser  power  output  has  been  shown  to  be 
correlated  to  the  presence  of  CN  species  in  the  plasma. 
The  chemistry  of  CN  in  the  CO  laser  will  be  discussed 
using  our  detailed  numerical  model  together  with  ex- 
perimental evidence  from  our  laboratories  and  else- 
where. As  a prelude  to  the  N2  chemistry  model,  the 
earlier  results  on  O2  additions  will  be  summarized. 

(l)W.L.  Morgan  and  E.R.  Fisher,  manuscript  in  pre- 
paration . 
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Experimental  measurement  of  the  1st  and  2nd  overtone  spontaneous  emis- 
sions spectrum  in  a low  pressure  CO/He  plasma  has  enabled  determination 
of  the  vibrational  distribution  up  to  v = 30.  The  variation  in  the  CO 
distribution  as  a function  of  CO  and  He  partial  pressure  and  plasma  current 
has  been  correlated  by  a detailed  model  of  the  steady  state  energy  transfer 
dynamics.  The  agreement  between  calculated  and  measured  distributions  has  gen- 
erated a consistent  set  of  rate  coefficients  for  CO/He  mixtures  for  a tempera- 
ture around  150  K. 
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Introduction 


* 


The  CO  iaser,  due  to  its  high  efficiency,  has  undergone  detailed  investi- 
gation since  its  first  demonstrated  performance  by  Patel  [1],  A great 
deal  of  effort  has  gone  into  understanding  and  modifying  the  operating  condi- 
tions to  both  improve  efficiency  and  enhance  certain  spectral  operating  regions. 
These  effects  are  often  achieved  by  altering  the  plasma  environment  in  which 
the  excitation  process  is  produced  and  by  altering  the  chemical  kinetic  pro- 
cesses affecting  excitation.  In  order  to  aid  our  understanding  and  analysis 
we  have  underway  a detailed  study  of  plasma  effects  which  are  observed  by 
analyzing  the  overtone  spectra  from  a flowing,  continuous-electric-discharge, 

CO  laser  system. 

The  1st  and  2nd  overtone  spectra  are  measured  using  the  end-light  spon- 
taneous emission  from  the  discharge  tube.  In  this  configuration  we  have 
been  able  to  observe  output  corresponding  up  to  the  35th  vibrational  level  in 
CO  with  substantial  signal-to-noise  [2].  We  are  not  yet  limited  in  observation 
by  any  physical  considerations.  Using  this  overtone  spectra  we  have  calculated 
the  CO  vibrational  distribution  up  to  V^$30  and  thereby  determined  directly 
the  effects  on  the  vibrational  distribution  resulting  from  changes  in  the 
plasma  environment,  either  electrical  and/or  chemical.  The  experiments 
analyzed  and  reported  in  this  paper  are  restricted  to  the  binary  mixture  CO/He. 
We  have  varied  the  relative  partial  pressures  of  CO  and  He  in  a systematic 
fashion  and  have  found  a set  of  energy  transfer  rate  coefficients  which 
are  both  consistent  with  available  rate  data  and  will  model  all  observations 
successfully.  This  paper  reports  these  rate  coefficients  for  the  CO/He 
system,  the  measurement  technique  used  to  generate  the  CO  vibrational  dis- 
tributions, and  the  analytical  laser  model  developed  to  correlate  the  observed 
distributions.  Measurements  coordinated  with  analytical  modeling  are  underway 
on  other  relevant  mixtures  involving  CO,  e.g.  CO/^/He,  C0/0,,/He  and  CO/Ar. 

The  results  of  these  studies  will  be  detailed  in  forthcoming  papers. 
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II.  Measurements 


The  overtone  spectra  in  CO  were  observed  using  the  spontaneous  end-light 
emission  from  a liquid-nitrogen  cooled  electric-discharge  laser  tube,  with  all 
| laser  action  blocked.  The  end  light  was  analyzed  by  a quarter-meter  low 

resolution  scanning  spectrometer  with  a cooled  PbS  detector.  The  data 
was  simultaneously  displayed  on  an  X-Y  plotter  and  recorded  on  an  FM 
recorder  for  computer  analysis. 

The  medium  was  checked  for  optical  thinness  by  placing  a mirror  at  the 
end  of  the  discharge  tube  to  reflect  the  radiation  back  to  the  detector.  It 
was  observed  that  the  increase  in  intensity  was  a constant  proportion  of  the 
signal  indicating  no  significant  self-absorption. 

One  of  the  features  of  this  tube,  which  later  played  an  important  role 
in  the  analysis,  was  a teflon  collar  placed  over  the  edge  of  the  anode  to  block 
out  plasma  contact  at  the  edge.  This  collar  insured  that  the  discharge  connected 
to  the  large  area  of  the  inner  surface  of  the  concentric  ring  electrode  as 
shown  in  Figure  1. 

The  gases  were  premixed  and  pumped  into  the  discharge  volume.  The  compo- 
sition could  be  accurately  controlled  using  a calibrated  rotometer  for  each 
constituent.  Starting  from  the  nominal  mixture  that  gave  optimal  laser 
performance,  the  1st  and  2nd  overtone  CO  spectra  were  measured  in  mixtures 
for  which  0.5  < P <1.5  and  4 < < 16  torr.  For  each  run  the  voltage 

across  the  positive  column  was  monitored  so  that,  together  with  the  composi- 
tion and  total  pressure  data,  the  plasma  could  be  characterized. 

The  intensity  versus  wavelength  data  was  transferred  to  a high-speed 
digital  computer  and  analyzed.  As  a first  step,  the  spectrum  was  corrected 
for  system  response  functions  by  using  the  recorded  spectrum  of  a calibrated 
"grey-body”  viewed  from  behind  the  end  of  the  laser  tube.  The  data  was 
simulated  using  standard  infrared  vibrational  transition  profiles,  convoluted 


with  the  spectrometer  window  profile  and  then  fitted  by  a least-squares 
adjustment  technique  to  the  observations.  As  a result  of  the  fit  we  obtain 
the  relative  population  distribution  in  CO.  Because  of  the  overlap  of  the 
high  v states  of  the  2nd  overtone  with  the  low  v states  of  the  1st  overtone 
spectra  the  least-squares  fit  had  to  be  done  iteratively. 

A representative  spectra  and  the  resulting  fit  are  shown  in  Figure  2. 

In  the  fitting  procedure  the  'goodness'  of  the  fit  was  affected  by  changing 
the  kinetic  temperature.  Although  this  had  little  affect  upon  the  relative 
population  distribution  it  did  affect  the  overall  quality  of  the  fit  to  a 
degree  that  permitted  us  to  obtain  the  rotational/kinetic  temperature 
to  within  + 15K.  The  resulting  temperature  could  be  checked  against  the 
temperature  predicted  according  to  the  profile  of  the  2-0  P-branch  contour 
which,  under  conditions  of  no  overlap  of  the  2nd  overtone  emission,  could  be 
obtained  unambiguously.  The  results  from  the  two  methods  are  in  agreement 
to  within  the  data  spread.  It  should  also  be  noted  that  in  our  experimental 
technique  the  minimum  signal-to-noise  in  the  first  overtone  spectra  occurred 
in  the  highest  vibrational  level  spectra  monitored,  35-33,  and  was  typically 
25  to  1. 

Ill . Steady  State  Laser  Model 

In  order  to  interpret  the  experimentally  determined  vibrational  distribu- 
tions, we  have  developed  a detailed  kinetic  model  incorporating  a non-Boltzmann 
electric  energy  distribution  code  [3]  and  including  the  dynamics  of  vibra- 
tional energy  transfer  in  a multi-level  system.  The  laser  code  uses  a 
rate  equation  formulation  and  includes  all  Av  = + 1 processes  for  vibration- 
translation  (VT)  and  vibration-vibration  (VV)  energy  transfer,  spontaneous 
fundamental  and  overtone  emission  processes  and  multi-quantum  electron 
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Figure  2:  Experimentally  observed  first  and  second  overtone 

spectrum  for  CO  from  a CO/He  mixture  (1/7)  at  8.5  Torr  total 
pressure  and  155  K temperature.  The  points  indicated  are  from 
the  computer  fit  to  the  observed  spectrum. 




excitation  and  de-excitation  (super-elastic)  processes.  The  typical  laser 
code  calculation  presented  here  includes  40  levels  in  CO  and  He  as  a diluent. 

All  electron  processes  are  generated  from  the  non-Boltzmann  electron 
energy  distribution  code  described  extensively  elsewhere  [3].  The  code  solves 
the  steady  state  Boltzmann  equation  using  a defined  E/N  and  effective 
vibrational  temperature  for  each  vibrational  level.  The  electron  code  is 
directly  coupled  to  the  kinetic  model  so  that  if  the  vibrational  temperature 
of  CO  changes  enough  to  influence  the  electron  energy  distribution  through 
super-elastic  processes,  the  electron  code  is  recalled  and  the  rates  recal- 
culated. Although  super-elastic  feedback  has  profound  effects  on  the  high 
energy  tail  of  the  electron  energy  distribution  [3],  the  relatively  low-energy 
electrons  responsible  for  vibrational  excitation  are  essentially  unaffected 
by  changing  vibrational  temperature.  Nonetheless  the  calculations  presented 
in  the  following  section  include  a fully  converged  super-elastic  contribu- 
tion to  the  final  CO  vibrational  distribution. 

The  rate  equations  describing  the  vibrational  populations,  in  the  steady 
state  approximation,  lead  to  a set  of  coupled,  nonlinear  algebraic  equations. 
These  equations  are  numerically  solved  using  a double  precision  Gauss  Elimination 
method.  The  laser  kinetic  model  and  the  non-Boltzmann  electron  energy  distribu- 
tion code  is  also  available  for  calculating  the  vibrational  distributions  for 
2 coupled  diatomics  with  arbitrary  diluents,  with  or  without  stimulated  emission 
processes.  Examples  of  the  latter  expanded  code  are  being  used  to  correlate 
C0/N_/He  experiments. 

Critical  to  the  accurate  prediction  of  the  CO  vibrational  distribution 
are  reliable  rate  coefficients  and  cross  sections.  The  cross  sections  for 


vibrational  excitation  by  electrons  used  in  the  electron  code  are  due  to 


W 

f 


Ehrhardt,  et  al  [4],  The  other  electron  cross  sections  incorporated  into 
electron  calculation  are  discussed  eLsewhere  [3]. 

The  A-coef ficients  for  spontaneous  emission  in  both  fundamental  and 
overtone  bands  were  obtained  from  earlier  work  [2]. 

The  VT  and  VV  rate  coefficients  used  in  these  calculations  are  an  opti- 
mized set  of  coefficients  based  on  available  experimental  rates,  theoretical 
predictions  and  analytical  modeling.  The  basic  features  of  the  VT  rate 
coefficients  follows  closely  the  calculations  of  Verter  and  Rabitz  [6] 
modified  by  v ' where  v is  the  upper  state  quantum  number.  These  VT 
coefficients  are  rather  similar  to  those  used  by  Rich  [7]  for  levels  1 < v < 20. 
Above  about  v = 20  the  VT  rates  used  here  show  a much  weaker  dependence  on 
increasing  v than  the  Rich  rates  [5]. 

The  VV  rate  coefficients  were  generated  using  the  sum  of  long-range 
and  short-range  potential  models  [8].  Adjustment  of  the  model  parameters 
was  performed  until  the  rate  coefficients  gave  reasonable  agreement  with 
available  experimental  rates  in  CO  as  shown  in  Figure  3.  The  scaling  of 
the  short  range  potential  model  to  upper  vibrational  states  was  accomplished 
using  Morse  oscillator  matrix  elements. 

j IV.  Results 

As  will  be  detailed  elsewhert  [51,  a consistent  set  of  energy  transfer 
coefficients,  i.e.  VT  and  VV,  was  determined  for  the  CO/He  system  in  the  temp- 
erature range  125  < T < 175K.  In  addition,  an  electron  density  and  E/N  were  chosen 
which  brought  the  predicted  stea  !y  state  CO  vibrational  distributions  into  ex- 
cellent agreement  with  the  measured  distributions.  A test  on  the  consistency  of 
the  choice  for  the  rates  and  the  plasma  parameters  should  consist  of  equally  good 
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VIBRATIONAL  LEVEL 


Figure  3:  Calculated  probabilities  for  W pumping  in  CO  at 

135  and  173  K in  comparison  to  available  experimental  data. 
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comparisons  between  predicted  and  measured  distributions  over  a wide  range 
of  CO/He  partial  pressures  and  plasma  currents.  Since  the  plasma  neutral 
temperature  is  determined  from  the  experimental  data  and  the  CO  and  He 
partial  pressures  are  directly  measured,  a model  calculation  can  be  made  if 
the  E/N  and  e-density  are  known.  The  voltage  across  the  positive  column 
and  the  total  pressure  are  both  measured  so  the  E/N  for  the  i*^  experiment 
can  be  scaled  from  the  "baseline"  experiment  by 


CD 


where  the  subscript,  b,  refers  to  the  "baseline"  parameters. 

The  electron  density  is  related  to  the  plasma  current  through 
R 2* 

I = f f V [e]  r db  dr 
J0  *0  ° 


(2) 


where  = electron  drift  velocity  and  [e]  = electron  density.  We  can  signifi- 
cantly simplify  this  relationship  by  assuming  the  drift  velocity  is  independent 
of  radius.  Therefore,  using  the  E/N  as  defined  in  equation  (1)  together  with 
the  gas  partial  pressures,  the  drift  velocity  can  be  calculated  from  the 
electron  energy  distribution  code  [3].  If  the  integral  is  then  replaced  by 
"effective"  values  we  find 


1 - ! V [e]eff-  Deff 


(3) 


where  [e]  is  assumed  to  be  the  electron  density  responsible  for  the  laser 
excitation  and  D ^ is  the  effective  diameter  of  the  laser  plasma  responsible 
for  transporting  the  current.  If  D ^ is  assumed  constant,  the  above  relationship 
would  provide  a basis  for  scaling  e-density  between  experiments. 

Figures  (4)  through  (6)  exemplify  the  experimentally  determined  vibrational 
distributions  in  comparison  to  the  calculated  distributions  based  on  constant 
D ^ (dashed  lines).  The  general  agreement  is  poor.  However,  when  the  electron 
density  is  left  arbitrary,  the  solid  lines  ir.  the  above  figures  result. 
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Figure  4:  Steady  state  CO  vibrational  distribution  for  a 

CO/He  mixture  of  (1/7)  at  155  K.  The  plasma  parameters 
used  were  electron  density  = 8.5(9)  and  E/N  = 1.5 (-16). 
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RELATIVE  CONCENTRATION 


° EXPERIMENTAL 

CALCULATED 
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Figure  5:  Steady  state  CO  vibrational  distribution  for  a 

CO/He  mixture  of  (l/4)at  155  K.  The  plasma  parameters 
used  were  electron  density  = 5.8(9)  and  E/N  = 2. 7 (-16). 
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Figure  6:  Steady  state  CO  vibrational  distribution  in  a 
CO/Ite  mixture  of  (1/16)  at  135  K.  The  plasma  parameters 
used  were  electron  density  = 8.9(9)  and  E/N  s:  1.4(-16). 


The  agreement  in  all  cases  except  Figure  (>  Is  excellent.  Figure  6 represents 
the  experiment  with  the  largest  He  partial  pressure  and  shows  a substantially 


larger  measured  population  in  the  high  vibrational  states  than  predicted  from  the 
model.  It  is  precisely  this  experimental  condition,  i.e.  large  He  pressure, 
that  should  correspond  to  maximum  cooling  and  lowest  gas  temperature.  The 
solid  line  is  the  result  of  a model  calculation  for  a 135  K temperature, 

20  K below  the  other  predictions.  The  agreement  at  this  lower  temperature 
is  found  to  be  excellent. 

Four  of  the  above  runs  were  made  under  identical  partial  pressures  of  CO 

and  He  but  with  increasing  plasma  current  from  5 to  15  ma.  If  the  predicted 

plasma  parameters  from  these  Figures  are  used  together  with  equation  (3)  to 

calculate  D c , Table  I results.  Also  shown  in  Table  I are  effective  plasma 
et  i r 

diameters  obtained  photographically  using  Polaroid  film.  These  optically 
derived  effective  diameters  are  in  very  close  agreement  to  those  obtained  from 
the  modeling  calculations. 

In  order  to  show  the  el  feet  of  using  previously  accepted  VT  rates 
for  CO/He,  a run  identical  to  RUN  1 (Figure  4),  above,  was  made  using  the 
VT  rates  of  Rich  [ 7 i . The  results  for  this  comparison  is  shown  in  Figure  7. 
Above  about  v=25  where  the  VT  rates  of  CO/He  control  the  shape  of  the 
distribution  the  deviation  is  seen  to  be  large  and  in  disagreement  with  the 
measured  values  in  these  levels. 
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Figure  7:  Comparison  of  steady  state  CO  vibrational 

distributions  calculated  using  the  VT  rates  in  this 
paper(5)  and  the  VT  rates  of  Rich(7)  from  the  SSH 
theory.  The  parameters  used  are  the  same  as  in  Figure  4 


Table  1 


Effective  Plasma  Diameter  in  a CO/He  Laser* 


I(  mA) 

[e]eff (cm  3) 

VD(cm/sec) 

obs 

Dgjf (cm) 

5 

85 

5.0*  + 

1.0 

1.0 

7 

li0 

4.56 

1.2 

1.1 

10 

1.310 

3.76 

1.5 

1.3 

15 

6.59 

5.16 

1.8 

1.9 

* [CO]  = 1 torr,  [He]  = 7 torr,  D = 2.5  cm. 
1 ‘ tube 

4 normalized  to  the  5mA  observation 
t 5.06  = 5.0  x 106 
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Discussion 


The  ability  to  experimentally  determine  and  analytically  correlate 
CO  vibrational  distributions  to  high  vibrational  states  demonstrates  an 
important  new  tool  to  investigate  the  energy  transfer  dynamics  of  high 
vibrational  states.  In  the  experiments  reported  here,  the  agreement 
between  experiment  and  predicted  vibrational  distributions  up  through 
v = 30  has  generated  a consistent  set  of  energy  transfer  rates  for 
CO/He  mixtures  in  the  temperature  range  near  150  K.  These  rate 
coefficients  have  been  shown  to  differ  significantly  from  past  values 
in  vibrational  levels  above  v = 25.  Further,  the  quantification  of 
strong  radial  gradients  in  the  electron  density  provide  important 
information  for  use  in  radial  diffusion  models  to  gain  greater  insight 
into  ambipolar  diffusion  and  wall  loss  processes  in  low  density  CW  laser 
systems.  Additional  studies  involving  other  inert  diluents  ( to 
investigate  VT  mass  scaling  effects  ),  other  wall  temperatures  ( to 
investigate  the  temperature  dependence  of  W and  VT  rates  ) and  additive 

(9) 

species  such  as  ^ ( to  study  ^-CO  W pumping  ) and  0o  (to  observe 
chemical  effects  in  hgih  states  ) are  underway. 
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A coefficients  for  spontaneous  emission  in  CO 

Allan  J.  Lightman 
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Edward  R.  Fisher 

Research  Institute  for  Engineering  Sciences.  Wayne  State  University.  Detroit.  Michigan  4H202 
(Received  6 July  1976) 

Simultaneous  end  light  emission  measurements  of  the  first-  and  second-overtone  hands  of  CO  have  been 
obtained  in  a flowing  CO-N,-He  laser  plasma  Interpretation  of  the  band  structure  using  A coefficients 
from  Young  and  Eachus  show  consistent  agreement  with  observation  up  to  V 26.  This  agreement 
supports  the  use  of  these  A coefficients  in  laser  and  kinetic  applications 

PACS  numbers:  42  60  Cz,  52.25  Ps 


Fundamental  to  studies  on  improved  CO  laser  sys- 
tems and  to  studies  using  the  CO  laser  as  a probe  for 
reacting  chemical  systems  is  an  accurate  knowledge 
of  the  Einstein  A coefficients  for  the  CO  molecule.  Since 
the  CO  molecular  laser  operates  on  fundamental  transi- 
tions up  to  very  high  vibrational  levels,  a knowledge  of 
the  scaling  of  these  A coefficients  with  increasing  vi- 
brational level  is  particularly  important  in  predicting 
laser  performance  and  excited  CO  vibrational  popula- 
tions.1 Current  CO  laser  studies2-4  scale  the  A coeffi- 
cients with  vibrational  levels  using  the  matrix  elements 
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obtained  from  the  cubic  dipole  moment  function  of  Young 
and  Eachus.  Although  there  has  not  been  a serious 
question  on  the  validity  of  this  scaling  procedure,  we 
provide  experimental  evidence  to  support  the  use  of  the 
Young  and  Eachus  vibration-rotation  matrix  elements 
for  scaling  the  .4  coefficients  up  to  at  least  r - 26. 

As  part  of  a program  aimed  at  characterizing  the 
effects  of  additive  gases  on  the  performance  of  dis- 
charge-sustained CO  molecular  lasers,  we  routinely 
observe  spontaneous  end  light  emission  from  overtone 
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FIG . 1.  Experimentally  observed  infrared  emission  from  CO 
in  a low-pressure  [lowing  molecular  laser  under  liquid-nitro- 
gen cooling. 


bands  of  CO  as  a means  of  determining  the  CO  vibra- 
tional distributions  present  in  the  laser  system  The 
simultaneous  observation  of  both  the  first-  and  second- 
overtone  emission  in  CO  permits  a consistency  check 
on  the  vibrational  distributions  and  .4  coefficients  avail- 
able for  the  CO  molecule. 

By  using  standard  spectroscopic  techniques6-8  and 
the  A coefficients  of  Young  and  Eachus,  we  have  deter- 
mined the  rotational  temperature  and  vibrational  dis- 
tribution under  both  liquid -nitrogen  and  dry-ice  cooling 
from  the  first -overtone  band  structure  up  to  r 26. 

From  the  observed  second -overtone  band  structure, 
relative  A coefficients  corresponding  to  Ar  =3  transi- 
tions in  CO  have  been  determined.  These  experimental 
second -overtone  .4  coefficients  agree  very  well  with  the 
calculated  values  of  Young  and  Eachus  over  the  entire 
range  of  vibrational  transitions  from  r 3 to  v - 26. 

This  agreement  supports  the  extensive  use  of  these 
coefficients  in  interpreting  reaction  rate  data'-11  and  in 
CO  molecular  laser  calculations. 1-4  This  letter  out- 
lines the  experiments  and  methods  used  to  reach  this 
conclusion. 

The  infrared  spectra  was  collected  as  end  light  emis- 
sion from  a flowing  15-Torr  CO/N2/He  mixture 
(1:2:10)  in  a laser  tube  2.5  cm  in  diameter  and  80  cm 
long.  Typical  gas  velocities  of  1 m /sec  were  used  under 
both  liquid-nitrogen  or  dry-ice/alcohol  cooling.  The 
cw  discharge  was  supported  by  annular  electrodes  at 


in 

z 


a 

LxJ 

cr 


50 

40 

30 


: k, 

; %«, 


. NTFNS'TT  if.XP) 
- intensity  irm 


V U*  °ry,ce  co< 

. .jiii.N-, 


or  ...  . 

Zi  2.5  Z7 


UC.  D 

■A  COOL ! NG 


2 9 3.1  3 3 

WAVELENGTH  (MICRON) 


FIG.  2.  Analytical  fit  to  first -overtone  CO  emission  under 
both  dry-ice/alcohol  and  liquid-nitrogen  cooling. 
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typically  8 kV  and  10  mA.  A J--m  Jarrell-Ash  mono- 
chromator with  a 2.  5-g  blaze  grating  and  dry-ice 
alcohol -cooled  PbS  photodiode  located  50  cm  from  the 
discharge  tube  was  used  in  collecting  the  spectra. 

Signal  processing  was  carried  out  with  a chopper.  PAR 
amplifier,  and  either  an  v-v  recorder  or  a multichannel 
analog  tape  recorder  system.  Analog  to  digital  conver- 
sion was  performed  for  data  analysis.  The  wavelength 
dependence  of  the  monochromator,  detector,  and  re- 
cording instrumentation  was  calibrated  by  replacing  the 
end  light  emission  with  a 900K  blackbody  source  The 
field  of  observation  was  apertured  to  the  central  10 
mm  of  the  laser  tube.  For  fixed  mixture  ratios  of  the 
laser  gases,  the  infrared  signatures  were  reproducible 


FIG.  .4.  CO  vibrational  distributions  extracted  from  the  ob- 
served first -overtone  emi  ~siun  in  Fig.  2 under  liquid-nit  rogen 
( . TtoX  l:i.*>  K)  and  dry  ice  alcohol  < . 2ln  K.1  cooling 
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FIG.  4.  Experimentally  obtained  second -overtone  A coeffi- 
cients for  CO  in  comparison  with  the  calculations  of  Ref.  5. 


to  better  than  5%.  A typical  first  and  second  overtone 
spectra  obtained  under  liquid -nitrogen  cooling  is  shown 
in  Fig.  1.  Using  the  data  obtained  and  the  exact  values 
of  the  A coefficients, 5 the  medium  can  be  shown  to  be 
optically  thin. 

The  rotational  temperature  was  determined  from  the 
R branch  of  the  2 — 0 first -overtone  band  by  a least- 
squares  fit  to  a Boltzmann  distribution.  Using  this  rota- 
tional temperature  and  A coefficients  for  the  CO  first- 
overtone  system  as  determined  by  Young  and  Eachus,5 
see  Table  1,  the  vibrational  distribution  was  calculated 
by  least -squares  adjustment  of  a profile  generated  from 
the  spectroscopic  data  of  Yardley.2  Typical  spectral 
fits  and  calculated  vibrational  distributions  are  shown 
in  Figs.  2 and  3,  respectively.  These  distributions 
show  the  typical  non-Boltzmann  form  expected  in  a dis- 
charge-supported CO  system. 1,4,13 


These  calculated  vibrational  distributions  and  rota- 
tional temperatures  were  then  least -squares  fitted  to 
the  observed  second -overtone  spectra  and  A coefficients 
were  generated.  A coefficients  for  the  second -overtone 
transitions  in  CO  up  to  r ...  26  obtained  at  both  liquid- 
nitrogen  and  dry-ice/alcohol  temperatures  are  shown 
.n  Fig.  4 in  comparison  with  theoretical  values.  The 
agreement  over  the  entire  range  of  accessible  vibra- 
tional levels  is  good. 

Since  a consistent  interpretation  to  both  the  first-  and 
second -overtone  spectra  has  been  obtained  using  the  A 
coefficients  of  Young  and  Eachus,  we  feel  this  agree- 
ment supports  the  use  of  these  A coefficients  in  kinetic 
and  laser  applications. 
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ABSTRACT 

Direct  measurements  of  the  vibrational  distribution  to  v^30  obtained 
from  the  1st  overtone  emission  from  CO  has  been  made  in  the  presence  of 
varying  amounts  of  O2.  These  observed  distributions  show  enhanced  vibrational 
populations  for  small  additions  and  show  upper  level  deoopulation  for  larger 
additions.  Both  observations  support  a recent  plasma-kinetic  model  of  the 
effect  of  O2  in  low  pressure  CO  laser  plasmas  (1). 
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The  effect  of  0^  on  CO  laser  performance  has  been  extensively  re- 
ported.^ ^ Until  recent  detailed  model  calculations , ^the  effect  of  O2 
on  laser  performance  has  been  misunderstood.  In  this  note  we  present 
direct  IR  measurements  on  the  effect  of  O2  on  the  CO  vibrational  distri- 
bution which  supports  these  recent  model  calculations. 

02  in  small  additions  is  known  to  significantly  enhance  laser  output 
performance  while  in  larger  additions  the  laser  output  is  known  to  degrade. 
This  behavior  has  been  related,  by  model  calculations,  to  plasma-chemical 
processes. ^ In  small  additions,  (i.e.  - 10%  CO),  0^  replaces  polymer 
C0+ions  as  the  dominant  positive  ion.  Due  to  the  reduced  recombination 
rate  of  0^  , the  electron  density  and  hence  the  vibrational  excitation 
rate  of  CO  increases.  In  addition,  the  reduced  E/N  (electric  field/total 
gas  density)  noted  for  small  02  additions  leads  to  reduced  plasma  heating 
and  hence  lower  plasma  temperatures.  In  larger  additions,  sufficient  O2 
dissociation  by  direct  electron  impact  is  noted  to  produce  a reduction  in 
CO  vibrational  energy  by  C0+  0 vibrational-translational  (VT)  relaxation. 

This  VT  relaxation,  as  well  as  other  associated  plasma-chemical  processes!^ 
leads  to  increased  plasma  heating  which  further  reduces  laser  output  power 
and  efficiency.  These  interpretations  are  collaborated  by  the  measurements  and 
model  calculations  presented  in  this  note. 

For  varying  O2  additions,  the  IR  1st  overtone  spectra  of  CO  has  been 
measured  directly  in  a liquid  Ng  cooled  laser  plasma^ ’^The  spectra  are 
recorded,  digitized  and  numerically  processed  to  obtain  CO  vibrational 
distributions  up  to  v » 30.  The  measurements  and  deconvolution  procedure 
typically  leads  to  an  accuracy  of  less  than  ±20%  in  the  CO  vibrational  dis- 
tribution. The  magnitude  of  this  uncertainty  is  directly  related  to  water 
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absorption  bands  in  the  laboratory.  For  pure  CO/He  mixtures,  previous  work 

has  shown  quantitative  agreement  between  measured  distributions  and  steady 

( 8 1 

state  model  laser  calculationsv  'supporting  the  energy  transfer  and  electron 

excitation  rates  used  in  these  calculations. 

The  steady  state  laser  model  used  in  these  calculations  incorporates  a 

(1  9) 

non-Boltzmann  electron  energy  distribution  code  ’ 'and  includes  all  Av  = ±1 , 
VV  and  VT  energy  transfer  processes.  The  steady  state  algebraic  equations 
characterizing  the  vibrational  distribution  are  solved  using  a Gaussian 
elimination  procedure.  The  CO/He  VV  and  VT  energy  transfer  rate  coefficients 
at  temperatures  around  150K  (the  temperatures  encountered  in  this  study)  are 
an  optimized  set  which  have  been  shown  to  give  agreement  between  measured 
and  calculated  CO  vibrational  distributions  over  a wide  range  of  partial 
pressures  and  plasma  conditions^8- ^ 

In  Figure  1 we  show  available  data  and  calculations  on  vibrational 
relaxation  (VT)  of  CO  and  by  0 atoms.  Using  the  recent  unpublished 
measurements  of  Lewittes,  et.  al  as  presented  by  Kelley  and  Thommarson, 
we  estimate  that  the  relaxation  time  of  CO  + 0 collisions  at  150  K is 
about  8 usec-atm. 

Calculations  on  the  effect  of  0 atom  VT  relaxation  in  high  levels  of 

CO  requires  a mass-scaling  effect  with  increasing  vibrational  level  over 

(8) 

previously  derived  CO/He  VT  rate  coefficients.  In  providing  this  mass- 

(12) 

scaling  we  have  drawn  on  the  SSH  theory  as  developed  by  Lattr  'and  scaled 
the  CO/He  scaling  with  increasing  vibrational  level  by  the  ratio  of  reduced 
collision  mass  as  appears  in  the  exponent  of  the  adiabaticity  factor.  Thus, 
the  CO-O  scaling  is  a factor  of  2.8  larger  than  CO-He  at  v=20,  4.6  larger  at 
v=30  and  6.4  larger  at  v=39.  A further  quantification  of  the  effect  of  col- 
lision mass  on  VT  scaling  is  underway. 
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Figure  1:  Vibration-translation  relaxation  data  for  0 atoms 

on  ^ and  CO. 
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The  results  of  measured  CO  vibrational  distributions  without  0^  and 

with  0.1  torr  added  0^  are  shown  in  Figure  2.  The  solid  lines  represent 

the  predicted  distributions  for  the  stated  plasma  properties,  i.e.  electron 

density  and  E/N.  The  E/N  of  the  plasma  is  scaled  between  runs  based  on 

measured  positive  column  voltages  and  total  gas  pressure  while  the  electron 

density  is  adjusted  to  bring  calculated  and  measured  distributions  into 
(8 ) 

optimum  agreement.  ' 

A set  of  measured  vibrational  distributions  for  larger  02  additions 

up  to  1 torr  is  shown  in  Figure  3.  Clearly  for  0^  partial  pressures  above 

0.1  torr  the  upper  vibrational  levels  of  CO  become  preferentially  deactivated. 

This  behavior  is  characteristic  of  an  enhanced  VT  process  such  as  that  due  to 

CO  + 0 collisions.  Based  upon  this  observation  we  have  generated  theoretical 

1 5 

distributions  for  varying  0 atom  densities  up  to  5.2  /cc  as  shown  in  Figure  4. 
The  qualitative  behavior  shows  excellent  agreement.  By  assuming  appropriate  0 
atom  densities  from  Figure  4,  the  agreement  between  measured  and  predicted 
distributions  becomes  quantitative.  The  derived  0 atom  densities,  0?  additions 
and  calculated  dissociation  fractions  are  shown  in  Table  1.  We  should  indicate 
that  this  identical  behavior  has  been  observed  in  the  presence  of  N^,  suggesting 
that  C>2  effects  in  CO/N^/He  mixtures  might  be  of  similar  origin. 


Table  1 


Run 

0 Atoms  (microns) 

OoDensities  (microns) 

°l  Dissociation 

B3 

10 

— C 

250 

2.2 

B4 

24 

500 

2.6 

B5 

36 

1000 

1.8 

100 
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Figure  2 


Measured  and  calculated  CO  vibrational  distributions  for 
a 1.5/12.5  torr  CO/He  mixture  at  155  K with  and  without 
Op.  □ corresponds  to  no  02,  electron  density  = 1.0(10), 

E/N  = 1 . 68 (-16)  and  a current  of  10  ma;  O corresponds 
to  0.1  torr  02,  electron  density  = 2.0(11),  E/N  = 1.02(-16), 
and  a current  of  10  ma. 


RELATIVE  CONCENTRATION 


Figure  3 
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Measured  CO  vibrational  distributions  for  a 1.5/12.5  torr 
CO/He  mixture  at  155  K and  10  ma  for  various  0?  additions 
The  symbols  +/0/D/A  correspond  to  0/0.2570.5/1.0 
torr  added  0„. 
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Figure  4:  Calculated  CO  vibrational  distributions  for  a 1.5/12.5  torr 

CO/He  mixture  at  155  K,  electron  density  = 7.5(10)  and 
E/N  = 1 . 24 (-16)  for  various  0 atom  additions.  The  symbols 
0/0  /A/-f-/X  correspond  to  0.8/8.5/26/42/83  microns 
of  added  0 atoms. 
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Several  observations  emerge  from  the  agreement  between  these  measured 
and  calculated  CO  vibrational  distributions  in  the  presence  of  0 atoms.  The 
first  is  that  the  general  quantitative  agreement  is  very  good,  lending  supp- 
ort to  the  rate  coefficient  values  used  in  the  calculations.  The  second  is 
that  the  magnitude  of  the  dissociation  required  to  exDlain  these  results  is 

consistent  with  our  past  ion-chemistry  calculations  as  well  as  the  reported 

(13) 

measurements  of  others  Lastly,  we  point  out  that  small  additions  of 

0^  not  only  leads  to  a reduction  in  E/N,  as  observed  previously,  but  also  to 
a substantial  increase  in  electron  density  at  constant  current.  This  obser- 
vation suggests  that  the  plasma  undergoes  a significant  radial  contraction  for 
small  O2  additions,  since  the  drift  velocity  remains  essentially  constant. 
Although  large  contraction  effects  have  been  noted  in  other  plasmas^  ^ , we 
do  not  yet  understand  the  plasma-kinetic  origin  of  this  present  observation. 

The  complicating  factors  of  polymer  ions,  radial  plasma  temperature  gradients 
and  homogeneous  versus  heterogeneous  electron  loss  processes  do  not  appear  to 
lead  to  a simple  explanation  for  this  contraction. 

Lastly,  we  would  again  indicate  the  importance  of  plasma-chemical  processes 
in  molecular  laser  plasmas  and,  as  in  the  case  with  the  0^  additive,  these  pro- 
cesses can  have  a dominant  influence  on  laser  performance. 
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